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1. Introduction 

Biregular automorphisms of rational surfaces with positive topological entropy present a major interest in complex 
dynamics, but their construction remains a difficult problem of algebraic geometry. For an overview of this problem, 
we refer to |i7| and to the references therein. It is a paradoxal fact that these automorphisms, although hard to construct, 
can occur in holomorphic families of arbitrary large dimension (see [2|). Besides, the automorphism group of a given 
rational surface can carry many automorphisms of positive entropy; we refer to and IS for recent and interesting 
results on this topic. In this paper, pursuing our previous work [7], we concern ourselves with families of automorphisms 
of rational surfaces and related rigidity properties. 

Let X = {Xt)teB be a holomorphic family of complex manifolds parameterized by a smooth base B. The theory of 
deformations allows us to associate to X a number m(X) of generic parameters between and dimZ? that measures 
the variation of the complex structure of the generic fiber of the family X (see |7|). We say that a biholomorphism / 
of a complex manifold X is rigid if any deformation X of X on which the automorphism / can be extended satisfies 
m(X) =0. If / is an automorphism of X that can be extended to a deformation X, then the Kodaira-Spencer map at 
the central fiber commutes to the action /* of / on (X.TX); therefore the image of KS(X) is fixed by this action. In 
particular m(X) < dimker(/* — id); this inequality allows to give on the one hand a upper bound for m(X) and on the 
other hand a simple criterion to say whether {X,f) is rigid: 

Theorem A. — Let f be an automorphism ofX and let f* be the corresponding action of f on H' (X,TX). If X is a 
deformation ofX such that f lifts to an automorphism ofX, then 

m(X) <dimker(/*-Id). 

In particular, if 1 is not an eigenvalue off* then f is rigid. 

An automorphism / of X is infinitesimally rigid if / cannot be lifted to an automorphism of a non-trivial infinitesimal 
deformation of X, or equivalently if 1 is not an eigenvalue of /* acting on H' (X,TX). Remark that an automorphism 
can be rigid without being infinitesimally rigid, but the converse holds. The main problem in order to apply TheoremjA] 
in practical examples is that the computation of the action of Aut(X) on H'(X,TX), this is more delicate than the 
action of Aut(X) on the Picard group of X. The reason for this is that the first action is not defined for birational 
automorphisms, whereas the second is. Our purpose in this article is to compute this action in various examples. 

The first case we deal with is the case of rational surfaces carrying an anticanonical curve. This class of rational 
surfaces is particularly relevant for us because automorphisms on these surfaces are well understood and studied (see 
for instance lUTl and ||8j|). The most significant result we prove in this setting is the following: 

Theorem B. — Let X be a basic rational surface with K|- < endowed with an automorphism f, and assume that 
there exists an irreducible singular curve C on X which is linearly equivalent to — K^. Then f is not infinitesimally 
rigid. Besides, X has no nonzero holomorphic vector field. 

Unfortunately, Theorem |B] doesn't allow to decide if / is rigid or not, all it says is that the rigidity of / cannot be 
proved using an approximation at the first order. 

The second case we deal with is the case of rational Kummer surfaces. Let £ be the elliptic curve obtained by 
taking the quotient of the complex line C by the lattice A = Z[j] (resp. A = . The group GL(2; A) acts linearly on 
the complex plane and preserves the lattice A x A; therefore any element M of GL(2; A) induces an automorphism /m 
on £ X E that commutes with the automorphism i (resp. (|)) defined by l{x,y) = (uCjiy) (resp. (^{x,y) = (jxjy)). The 
automorphism lifts to an automorphism /m on the desingularization X of (£ x T,)/ {i) , which is a rational Kummer 
surface. Besides the topological entropy of /m is twice the logarithm of the spectral radius of M. 

Theorem C. — For any matrix M of infinite order in GL(2; Z[j]), the automorphism fM is rigid. 

The situation is drastically different for the lattice Z[i]. Indeed, we have the following statement: 

Theorem D. — For any niatrix M in GL(2;Z[i]) such that (detM)^ = 1, the automorphism fM is not infinitesimally 
rigid. However, there exist elements in GL(2;Z[i]) whose associated automorphism is infinitesimally rigid. 
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Note that the automorphisms of rational Kummer surfaces appear in different important problems. Let us for instance 
recall the following result ||6|: let / be a birational map of a projective surface that preserves a holomorphic foliation iT. 
If (deg/")„ grows exponentially, then up to birational conjugacy and to a finite covering, either / is a Anosov auto- 
morphism on C^/ (A X A) with A = Z[j] or A = Z[i], or / is a monomial map on x Moreover ^ is linear and / 
preserves at least two foliations. As we will see in ^3.2|these examples also appear in the so-called Poincare Problem. 



Lastly, if / is an automorphism of a basic rational surface X, we explain how to compute in practice the action of / 
on (X,TX) from the data of the sequence of blowups of P^ defining X and the explicit formula of / considered as 
a birational map of P^. For doing this, we develop specific theoretical tools, such as 1-exceptional divisors. Then we 
carry out completely the computation in one specific example, which is the one constructed in [7, Thm 3.5]; it is an 
automorphism obtained by blowing up 15 successive points of P^, which correspond to three infinitely near points of 
length five. Some of the computations are computer-assisted (with Maple), we provide as an appendix of this paper a 
commented Maple file. 

Let us describe in more details the organization of the paper. In ^ we deal with infinitesimal rigidity of rational sur- 
face automorphisms. We introduce in ^2.1 [ the notions of rigidity and infinitesimal rigidity, and we give in Theorem |2.2| 
a bound on the number of parameters of a family of rational surfaces endowed with an automorphism. This gives in 
particular a cohomological criterion to prove the rigidity of an automorphism. In §2.2 we study automorphisms on 



rational surfaces X whose anticanonical divisor is effective, and we compute in Theorems 2.4 and 2.5 the action of the 



automorphism group on (X,TX^ in specific situations. In particular, we prove that a rational surface automorphism 
is never infinitesimally rigid if there exists a singular irreducible anticanonical curve on the surface. In ^ we deal 
with automorphisms of rational Kummer surfaces. Section ^3.1 [ consists in recollections about the construction of these 
surfaces, as well as formulae for the action of the automorphisms on the Picard group. We limit ourselves to the case 



of hexagonal lattices. Finally, using the Atiyah-Bott fixed point formula, we compute in Corollary 1 3. 5 
the representation of Aut(X) in H' 
and 



(X,TX] when X is a Kummer surface of hexagonal type. In 5 3.2 



the character of 

following nni 

his approach is 



15 1, we present another approach of the construction of rational Kummer surfaces via foliations 
geometric in nature and allows to prove in an elegant way that these surfaces are indeed rational. We have tried to make 
a self-contained exposition of the theory for the reader's convenience, but we claim no originality in these two sections. 



In § 3.3 we compute explicitly the characteristic polynomial of the action of an automorphism of a Kummer surface on 



its infinitesimal deformation space: for hexagonal lattices, this is Theorem 3.7 which refines Corollary 3.5 for cubic 



lattices, this is Theorem ] 3. 9 1 The purpose of ^is the construction of a theoretical framework to compute effectively in 
any concrete example the action of a rational surface automorphism on the space of infinitesimal deformations of the 
surface. The main idea is to represent the space H' (X.TX^ of infinitesimal deformations of X by global sections of the 

TX(D) 

sheaf ' for some special divisors D called 1-exceptional divisors, they are studied in 54.1 Then we construct two 



TX 

different bases of 



TX(D) 
TX 



in 54.2 



and 5 4.3 



both will be useful for practical computations of the action of automorphisms. 
The last section is devoted to an explicit example, and relies on the machinery developed in ^ The strategy is explained 



in 5 5. 1 and the computations are sketched in 5 5.2 and 5 5.3 



Acknowledgments. — The first author would like to thank Dominique Cerveau and Alcides Lins-Neto for pointing out 
and explaining lITOllTSl . Both authors would like to thank Philippe Goutet for the pictures and for the nice presentation 
of the appendix. 



2. Infinitesimal rigidity 

2.1. Infinitesimal rigidity is a spectral property. — Let X be any complex manifold. An infinitesimal deformation 
of X is an equivalence class of deformation at order one of a holomorphic cocycle defining X. These deformations are 
classified by the cohomology group (X,TX) via the Kodaira-Spencer map. 

Recall that in IT] §5] we have associated to any local holomorphic deformation X a number of parameters m(X) that 
measures the variation of the complex structure of the generic fiber of X. 
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Definition 2.1. — An automorphism / of X is rigid if any deformation 3iofX over a smooth base such that / lifts to 
an automorphism of 3£ satisfies m(3£) = 0. 

If Z is a rational surface without holomorphic vector field, then X satisfies m(X) = if and only if X is locally 
holomorphically trivial (because X admits a universal deformation over a smooth base). 

Theorem 2.2. — Let f be an automorphism ofX and let f* be the corresponding action of f on (X,TX). IfX is a 
deformation ofX such that f lifts to an automorphism of3i, then 

m(3e)<dimker(/*-Id). 

In particular, if 1 is not an eigenvalue off* then f is rigid. 

Proof. — Let be a point in the base B of the deformation and let fi, be the automorphism of X lifting /. Since f/, is a 
morphism of deformations, the Kodaira-Spencer map KS/, : I^B — > H' (X^, TX^) satisfies the identity o KS/, = KS/,. 
Therefore the image of KS^ is contained in ker(f^ — Id). Let us recall that the rank of KS^ is equal to m(X) if b is 
generic. Thus, for generic b, we have the inequality m(X) < dimker(fjj — Id). The function b i— )• dimker(f| — Id) being 
lower semi-continuous, we get that the inequahty is vahd for all in fi. □ 

We put the following definition: 

Definition 2.3. — An automorphism / of X is infinitesimally rigid if / can not be lifted to an automorphism of a 
non-trivial infinitesimal deformation of X, or equivalently if 1 is not an eigenvalue of /* acting on 

Note that an automorphism can be rigid without being infinitesimally rigid, but the converse holds: if / is infinitesi- 
mally rigid, then it is rigid. 

2.2. Automorphisms of anticanonical basic rational surfaces. — Anticanonical surfaces are surfaces whose anti- 
canonical class is effective. Let X be a basic rational anticanonical surface, and let D be an effective divisor representing 
the class — Kx. Let us write X as a blowup 

X = Xn — > Xf]^\ — !■ . . . — > X\ — > Xo = P^, 

and let 7ti,...,7tAf be the blowup maps and Ei,...,Ef^ be the exceptional divisors. Then there exists a sequence 
(Z),)o<(<Ar of effective divisors such that each Di is an anticanonical divisor on X,, Djv = D, and for any integer i 
between 1 and A^, D, = 7t*A_i -Ei- 

An anticanonical curve on is a cubic. Then it is easy to prove that if D is reduced, then D is either an irreducible 
reduced curve of arithmetic genus 1, or a cycle of smooth rational curves, or one of the three exceptional configurations 
shown in the picture below: 




If D is reducible but non reduced, the situation gets more complicated, even if we know that the irreducible compo- 
nents are still smooth rational curves. We start with the case of irreducible divisors. 
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Theorem 2.4. — Let X be a basic rational surface with K|. < endowed with an automorphism f, and assume that 
there exists an irreducible curve C onX which is linearly equivalent to — Kx- Let Pf, Qf and 6c denote the characteristic 
polynomials off* acting on Pic{X), H' (X,TX) an(iH*'(C,N*^^) respectively. 

(a) IfCis singular, Qf = Pf x Qc and f is not infinite simally rigid. 

(b) IfC is smooth, let us denote by C,c the action off* on the complex line H^(C,i2^). 

Then Qf = Pf x Qc ^ ^ f ^'^^ infinitesimally rigid. 

In both cases, X has no nonzero holomorphic vector field. 
Proof. — Let us consider the long cohomology exact sequence associated with the short exact sequence 

— > iil(-C) — >^x — > ^x|c — ^ 

of sheaves on X. Since K|. < and since C is irreducible, | — Kx | = {C} so that / fixes globally the curve C. According 
to Serre duality for < 7 < 2, 

W{X,TXY :^U^-j{X,Cl],^Kx) {X,£lj, {-€)). 

On the other hand, by Hodge theory. Pic (X) (^zC ~ aiidH2-0(X) =H'''2(X) =0. Besides, h2(X,TX) =0 

andh^(Z,i2^) ~ h'^(Z,TX(g)Kx) = h^(X,TZ(— 3)) =0. Hence we have an exact sequence of /-equivariant complex 
vector spaces 

— ^H''(C,i2^|c) — >n^{X,TX)* — ^Pic(X)OzC^H'(C,n^|c) — >U^{X,TXy — >0. (2.1) 
Even if C is singular, we have a conormal exact sequence 

where the injectivity of the first arrow holds because C is a local complete intersection inX. Since the arithmetic genus of 
C is 1, the duaUzing sheaf COc is trivial. Besides, the conormal bundle N^^^ has degree — which is positive, so that by 

Riemann-RochandSerreduaUty,hO(C,N*/^) =-K| andhi(C,N^/^) =0. Therefore H^ (C,ni.|c) ~Hi(C,il^) ~C. 
It follows that 7t is the restriction morphism from H' to H' (C,i2^). 

Assume that C is singular. According to the adjunction formula, 

deg Q-c = deg ^x^c ~ ^c/x — ^x\c ~ '^^S ^c/x ~ ®c = 0. 

Besides, fi^ is non-trivial since C is singular Hence we obtain that h*' (C, 12^) vanishes (otherwise we would have an 
exact sequence — > Oc — > — ^ — ^ ^ where T is a torsion sheaf so that deg = deg T would be positive) 
so that by Riemann-Roch h^ (C, H^) vanishes as well. Therefore h' (C, ii^j^) vanishes, so that h" (Z , TX) = 0, and we 
obtain the /-equivariant exact sequence 

— )• h" (C, Nc/x) — > {X,TX) * — 5^ Pic(X) 0z C — 

which gives Qf = Pf xQc. The canonical class in H' (X,Q,x) is invariant, 1 is thus always an eigenvalue of Pf, and an 
eigenvalue of 2/ as well. 

Assume now that the curve C is smooth. Then H' (C,i2^) is canonically isomorphic to C. Besides, the restriction 
morphism 7t: H^ {X,Q.l^) — > (C,n^) ~ C is surjective since Jt([Kx]) = K^.C = -C^ > 0. Therefore h^{X,TX) = 
and we have two /-equivariant exact sequences 

— >}l°{C,Q}xic) — >ll\X,TXy — )-Pic(X)(8)zC — >C — >0 
0^hO(C,N*/x) ^H0(C,i2i|c) ^U%C,£lh)^0. 

This yields 

e/x(x-l)=/'/xecx(x-Cc). 
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Therefore, is an eigenvalue of Q /. Since belongs to the finite set { 1 , — 1 , i, — i, j , } , we have Z,^ = I. Hence 1 
is an eigenvalue of 12 . □ 

We now investigate the case of anticanonical divisors given by a cycle of smooth rational curves and leave the three 
other exceptional configurations to the reader 

Remark that if / is an automorphism of a rational surface with positive entropy, then | — Kx | is either empty or 
consists of a single divisor (otherwise / would preserve a rational fibration). In this last case, by replacing / by an 
iterate, we can assume that all irreducible components of the divisor are globally invariant by /. 

Theorem 2.5. — Let X be a basic rational surface without nonzero holomorphic vector field and assume that there 
exists a reducible and reduced effective anticanonical cycle D = Y!i=\Di on X such that Df < 0/or all i. For any 
automorphism f ofX leaving each Di globally invariant, let Pf, Qf and 9, denote the characteristic polynomials of f* 
acting on Pic{X), ll^(X,TX) ami H''(D;,N^yj^(— 5,)) respectively, where Si = sing{D)r\Di. Then 

* {x-iy 

Proof. — Let S be the singular locus of D. Then we have an exact sequence 

Since A^^y;^ has positive degree, H^(D,-,n^|^ ) ~ ~ C. Besides, since H''(D,-,n]).) = 0, we obtain the 

isomorphism H° (a, £^x|£,.) =i H''(A,A^* ./^)- It follows that 

1=1 

and that the map from ©[^j H''(D;,nJ[.|^ ) to H°(5,n^|5) is onto. Therefore 

1=1 (=1 
We now write down the /-equivariant exact sequence \2.\) used in the proof of Theorem |2.4[ 

r 

— ^ H''(D,A?*./;^(-5,)) — ^ H' (X,TX)* ^ Pic(X) ®z C — ^ C H''(X,TX)* 0. 

(=1 

Since we assumed that there were no nonzero holomorphic vector fields on X, we get Qf^- (x— 1)' = Pf x XYi=\ 9f Q 



3. A geometric example: Kummer surfaces 

We start by fixing some conventions and notations. For any group G and any vector space W, a G-module means 
a contravariant representation of G in W. For any element g in a group G, we write {g) for the cyclic group in G 
generated by the element g. We denote the vector space by V . If nothing else is specified, we consider V as a natural 
GL(2; A)-module, where a matrix M acts by M^^ Except explicit mention of the contrary, the letter A always denotes 
the hexagonal lattice Z[j] in C^, where j = e^'"/-'. 

3.1. Construction, Picard group and fixed points. — Let £ be the elliptic curve obtained by taking the quotient 
of the complex line C by the hexagonal lattice A and let A be the abelian surface T, xT,. The group GL(2; A) acts 
linearly on and preserves the lattice A x A. Therefore any element M of GL(2; A) induces an automorphism on A 
that commutes with the automorphism (|) of A of order three defined by — (jjc, jy). It follows that induces an 

automorphism of the singular surface A/ ((])) . After blowing up the singular points, we obtain an automorphism /m on a 
smooth model X, which is a rational Kummer surface. 
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Let us remark that (|) and (|) have the same 9 fixed points (we will see that these points also appear in § 3.2 1: 



/?i = (0,0),/?2= 0, 



P6 = 



3' 3 3 



,P3 



0,: 



1 



1 , 2j 



,/^8 



,P4 



01P5 



j 2 

3'3 



1 , 2j 2 



3 ' 3 ' 3 ' 3 / ' ^^'^ V 3 



1 2j 1 2j 

3 '3 3 



In particular the group G = ((])) ~ Z3 has 9 fixed points on A. We denote this finite set by S, it is a subgroup of the 
torsion points of A of order 3. Note that the group GL(2; A) leaves the set S globally invariant. We denote by H the 
kernel of the natural group morphism from SL(2;A) to SL(2;F9) where F9 = Z^lj]. This is a subgroup of finite index 
in GL(2; A) which fixes every point in S. 

There is a slightly different way to construct X: first we blow up the set 5 in A = E x £ and denote by A the resulting 
surface and by Ej the exceptional divisors corresponding to the points p,. The group G acts on A, and the quotient A/G 
is exactly our Kummer surface X. Then we have the following diagram, where 5 is the blowup map and K denotes the 
projection map from AtoX: 




5 

A X 

This diagram is equivariant for both actions of G and GL(2; A). We set £, = 7i(£',). Then a straightforward computation 
shows that K is the projection of the cyclic cover of X of order 3 along each rational curve Zs,-. 

Let us recall thatX is obtained by blowing up the projective plane in 12 distinct points. The number of blown up 
points can be determined easily by computing the Euler characteristic of X: indeed, x{A) ~ x{BlsA) = x{A) +9 = 9. 
As #G = 3, we get x(A \ {£,},) = 3x{X \ {£,},) . But x(A \ {£,},) = -9 so x{X \ {£,},) = -3 and 

X(X)--3 + 18 = 15=x(P^) + 12. 



Another way to see it is given in 5 3.2 



The Picard group of X can be described in the following way: since X is rational, Pic{X) ~ H^{X,Zx)- First we 
compute: 

h2(A,Z^) ~ h2(A,Za) ® (^0Z[£,] j ^ a|(A* X A*) ^0Z[£,-] j , 
where A* denotes the dual lattice of A. Then we get an isomorphism of GL(2; A)-modules: 

9 _ 

Pic(X) n\A,Z^f ^ (a|(A* X a*))"" ® Z[£,-]. (3.1) 

1=1 

It is easy to see that (a|(A* x A*))*^ is a free Z-module of rank 4, so that Pic{X) is a free Z-module of rank 13 (which 
is the number of blows up plus 1, as expected). 

LetM be GL(2;A) with eigenvalues a and p. Then the complex eigenvalues of the endomorphism A^M of the real 
vector space A^V are 

|a|2, |p|2, ap, ap, ap and ap. 

A straightforward but lengthy computation shows that the eigenvalues of A^M on the subspace ( Ag V)*^ are exactly 
|ap, IPp, ap and ap. Thus, for any element M in GL(2;A) with spectral radius r^, the spectral radius of M acting 
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on (a|(A* X A*))*^ is r^. This means that the first dynamical degreej^ of /m is given by the formula A,i(/m) = r^- 
More precisely, the characteristic polynomial of acting on Pic{X) is 

(x-l)^(x-|ap)(x-|p|2)(x~ap)(x~ap). 

In the remaining part of this section, we compute the number of fixed points of /m for any M in H. As a consequence, 
we compute the trace of the endomorphism with the help of Atiyah-Bott fixed point formula. We will obtain a much 
finer result in the next section (Theorem |3.7|i. 



In the sequel, we divide the set of fixed points of /m into two parts: the first part 0i consists of fixed points outside 
the exceptional divisors Ei ; and the second part ©2 consists of the remaining fixed points. 

Proposition 3.1. — Let M be a matrix in H such that ^ Id. Then the automorphism /m has |Tr(M) p + 11 fixed 
points on the Kummer surface X. 

Proof. — The first step consists in proving that the fixed points of /m are non-degenerate, i.e. that 1 is never an 
eigenvalue of the differential of /m at a fixed point. We deal with fixed points in ©i and ©2 separately. 

The points of ©1 correspond to points p in A\S such that fuip) lies in the orbit G ■ p, modulo the action of G. For 
any such point, the differential of /m identifies with M, jM or j^M, so that 1 is never an eigenvalue. 

The set ©2 can be described as follows: the automorphism /m acts by the projective transformation P(M) on the 
curve Ei via the identification 

£,— p(Tp,y)~P(y). 

Therefore we get two fixed points qa and on each Ej corresponding to the eigenspaces of M. This gives a concrete 

description of ©2, which consists of 18 points (two distinct points on each £,). To compute the differential of /m at such 
a fixed point, we can assume without loss of generality that this point lies in Ei. Let (e,/) be an eigenvector of M for 
the eigenvalue a and assume for simplicity that e 7^ 0. If {x,y) are the standard coordinates on V, we define coordinates 
(m, v) in A near Ei by putting x = m and y — uv. Then if we set Z = u^ and T = v, (Z, T) are holomorphic coordinates 

on X near q^. In these coordinates, q^ = ( 0, — ) and is given by 



/ , ni c + dT 

{ZJ)^[z{a + bT)\-^^ 

so that the eigenvalues of d/^ at qa are and a^^, which are different from 1. 

We can now apply the Lefschetz fixed point formula. As M) = — 0, we get that the number of fixed 

1: the trace on the part (a|(A* x A*))*^ is the 



points of fM is equal to Tr/^jj^jj^fj^j +2. Then we can use formula (|3.1 
sum 

|a|2 + |p|2 + ap + ap = |a + p|2; 

and since M belongs to H, /m acts trivially on the factor ©^=1 ^[Ej] so that the trace on this factor is 9. This gives the 
required result. □ 

We will now study in greater details the set ©1. Let us define three subsets Si, S2, S3 of A as follows: 

Si^{peA-~.S\fMip)^p}, S2^{peA-~.S\fM{p)^m}, S3^{peA-~.S\fMip)=<^Hp)}- 

To compute the cardinality of the 5,'s, we use the following statement: 

Lemma 3.2. — For any matrix P in GL(2; A) such that I is not an eigenvalue ofP, the automorphism fp of the complex 

I 1 2 

torus A has 1 — Tr (P) + detP fixed points. 



1. Let / be a bimeromorphic map on a Kahler surface X. The first dynamical degree of / is given byXi(/)= lim (|(/")*|)'/" where/* is the 
induced map fromH' '(X,R) into itself. 
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Proof. — Since the fixed points of fp are non-degenerate, we can apply the holomorphic Lefschetz fixed point for- 
mula 19] p. 426]; this gives 

1 



1 - Tr (P) + detP #Fix {fp) x 



l-Tr(P)+detP 



□ 



Corollary 3.3. — Let M be a matrix in H such that M ^ Id. One has the following equalities: 
#5i = |Tr(M)-2|^-9 and #^2 = #^3 = 1 1 +Tr (M) 1^ - 9. 



Proof. — This follows directly from Lemma [T2| since 

5i =Fix(/M)\5, 52 =Fix(/j2^) \5 and ^3 = Fix (/jm) \ 5. 

□ 

The group G acts transitively on each set Sj. Then ©1 can be written as the disjoint union of the quotients Si/G. 



Remark that we can obtain in this way the result of Proposition 3.1 indeed, if f = Tr (M) 



I i2 / I i2 

k-2 / i+H 



2 



#e. = ^-3 + 2|^^-3jH^I -7 
Since #©2 = 18, this gives #Fix(.^) = |f |^ + 11. 

As H° (X, TX) H2 (X, TX) = 0, the holomorphic Atiyah-Bott formula ffl Thm. 4. 12] yields the following result: 

~ Tr d/^ ' 

Proposition 3.4. — The trace of f^ acting on H' (X,TX) is equal to the sum — ^ d t(Id " df ) '^^^^^ through 

the fixed points of fM- 

Corollary 3.5. — Let M be a matrix in H such that ^ Id. Then the trace of f^ acting on H' (X,TX) is 8Tr(M). 
Proof. — We divide the fixed point set of fM in four parts: Si/G, S2/G, S$/G that form ©1 and ©2. We start by 



computing the contribution of ©2 in the sum of Proposition 3.4 Let us put t = Tr(M). Any pair of fixed points in a 
divisor £, yields the term 



a ^ + + a ^ 



(l-a3)(l-a-2) (l-a-3)(l-a2) 
6 3 

so that the contribution of 02 is —9t — - — - — j^p^- Now the contributions of S\/G, S2/G and S-^/G are respectively 



\ 3 / 2-f' \ 3 \+t \ 3 y 1+f 

3f 3f 

so that the contribution of ©i is f + ^ — - — j-j-y . Adding the contributions of ©i and ©2 we get — 8f. □ 
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3.2. An alternative construction via foliations. — In |Il0l[lll the authors consider the pencil fp4 — {!Fa}aeP^ 

where is given by the 1 -form co + ar| with 

&— [x^ ~ l)xdy— {y^ — l)ydx, r| = (jc^ — l)xdy— (y^ — l)x^ dx. 

The tangency set of fp4, given by co A r| = 0, is an algebraic curve; more precisely it is the union of 9 invariant lines 
that intersect at 12 points. Let C be the union of these points and lines. Let {x^y) be a system of coordinates of and 
let {Fa}a be the pencil of linear foliations on where Fa is described by Q-a = dy — adx. Let us denote hy p: — >A 
the canonical projection and by the pencil of foliations {(ja}a = {p*Fa}a- Let us consider the following elliptic 
curves on A: 

£o,i={0}x£, £i,o = £x{0}, "Eia = {{x,x)\x e T.} , T.i,j = {{x,-jx)\x e <£)} . 
Set = {2<),ij 'Ei,0, Eij}. Given F in'^ and m in Fix((|)) denote F,„ —F+m. The set 

£ = {F,„\m(EFix{^),Fe'^} 

consists in 12 elliptic curves Ei, . . ., £12; two such different elliptic curves intersect only at 3 fixed points of Fix((|)). 
Let C* be the union of Fix((|)) and these 12 elliptic curves. 

According to McQuillan |5, Chapter 8, Theorem 1] there exists a rational map A such that i|/*(fp4) = Q^. 

Let us explain the construction of Since the restriction of the morphism n of degree 3 to Ej is a biholomorphism, K 
maps any Ej to a rational curve of self-intersection —3. Besides K maps each elliptic curve K*'Ei, T.i e , to a rational 
curve with self-intersection —1. Denote by 7ti : X — > the blowing-down map of the curves £1, . . ., £12; one can 
see that Xq — {see LlSj Lemma 3]). The map V|/ is given by 

\|/ = 7ti7t5"': A--^pI 

Set (€h — Fix((|))* — \(/(Fix((|))) and let C* be the union of (f* and Fix((|))«. In lHOl Proposition 1] the author proves 

that C = C up to an automorphism of P^. In particular X is obtained by blowing up P^ in 12 points. 

The family {!J^a}a plays a very important role in the study of foliations; it allows to give a negative answer to the 
two famous following problems |10|: 

• Poincare Problem: "is it possible to decide if an algebraic differential equation in two variables is algebraically 
integrable (in the sense that it has a rational first integral) ?" llT2lll4lll3l 

• "Is it possible to recognize the genus of the general solution of an algebraic differential equation in two variables 
which has a rational first integral ?" 

Let M be an element of GL(2; A) with two eigenvalues fji, fj2 such that | < 1 < \/j2\- The linear transformation 
induced by M preserves the pencil {Fa}a and fixes ^^pi/y,, where (Pi,Yi) and {^2,12) ^"e the two eigenvectors 

of 'M; as a consequence /m preserves 7t*5*/?H.^|3,/Yi and K*d* p^,Fp^/y^. 

3.3. Action of the automorphism group on infinitesimal deformations. — In this section we study the GL(2; A)- 
module li^{X,TX) where X is the Kummer surface associated to an hexagonal lattice (we also give the results for 
cubical lattices at the end of the section). For this we will use two fundamental exact sequences: 

First exact sequence relating TA and TA via the blowup map 5 : 

_ 9 

O-^TA^5*TA^0t£,,T^.(-l)^O (3.2) 

!=1 

This sequence is GL(2; A)-equivariant and G-equivariant. Remark that T^j (— 1) is isomorphic to 0£, (1), but this iso- 
morphism is not canonical and cannot be made compatible in any way with the action of GL(2; A). 

Second exact sequence relating TA and TX via the cyclic cover 7t : 

_ 9 

{n^JAf ^ TX — > 0t£,,N£,/x (3.3) 

!=1 
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This sequence is GL(2; A)-equivariant. 

9 

Let us denote by the sheaf ^^^Ei*^E^^^)- ^'^^'^^ prove all technical results we need. 



Lemma 3.6. — The following assertions hold : 

(a) For any integer i with <i <2, the natural morphisms 

fH'(A,TA) — >H'(A,5*TA) 

\w\x,%J:A) — ^H'(A,7i:*7i,TA) — ^H'(A,TA) 

are Gh{2; A)-equivariant isomorphisms. 

(b) The cohomology groups (A, ^) ^ and (A, ^) ^ vanish. 

( c) The cohomology group H' (X, (71* TA)*^) vanishes. Besides, (X, (ti^TA)'^) is isomorphic to V as an H-module. 

(d) For \ < i <9, we have an isomorphism ofH-modules between H' (^EijNg./x) cmd V. 

Proof. — (a) We write the Leray spectral sequence for the pair (5*TA,A): we have 

E'^'" ^ H" (A, R''5h. 0^ (g) TA) and = Gr''H''+^ (A, 5*TA) . 

Since 5 is the projection of a point blowup, R^5*(>j — for ^ > so that the spectral sequence degenerates and we 
get H''(A,TA) = E^'" ~ E^,'" = H''(A,5*TA). The argument is the same for the second morphism: % being finite, 
R«Jt*TA = Ofor^>0. 

(b) The vanishing of H'(A,ir)'^ is straightforward as H^(A,i7') — 0. We see that H^{A,tJ^) is isomorphic as a 
G-module to the direct sum 0/LiTp,A, where (|) acts by the inverse of its differential at each point Therefore, 
H°(A, iT) ~ and the result follows since = {0}. 

(c) According to (a) and since G is finite, we have isomorphisms 

H' {X, {K,TA f) ~ H' {X, (7I*TA)) ^ ~ H' (A, TA) . 



Using (3.2i combined with (b) and (a), we obtain GL(2; A)-equivariant isomorphisms H'(A,TA)*^ ~ H'(A,TA)*^ for 



/ = 1, 2. Now H' (a, TA)'^ is isomorphic toy* (g)y (resp. A^V*(SV) for/ = 1 (resp. i = 2), where a matrix M in GL(2; A) 
acts by 'M(g) (resp. A^'M^M"'). This proves that H' (A, TA)"^ vanishes and that H2(A,TA)^ = H2(A,TA) 

since j^j^' = 1. This yields the result. 

(d)lfW — TpA, we have the Euler exact sequence on £, : 

^ 0£,(-l) ~^W®Oe, T£,(-l) 

which proves that the canonical bundle of Ei is canonically isomorphic to detlV* (g)c Oe, (^2). Since N^./^ is canonically 
isomorphic to C'£.(— 3), we get by Serre duality a chain of //-equivariant isomorphisms 

H' (£,-,N£,/x) ~ H^(£/, 0£,(3) (gK£,)* ~ detW«)H°(£,-, Oe,{1)Y detW(g)W:^ dety ®y 
This finishes the proof of the statement, since dety is trivial as an //-module. □ 

Theorem 3.7. — For any element M in H, the characteristic polynomial of the endomorphism f^ o/H' (X,TX^ is 
equal to (x^ — tT{M) x + l)^ . 

Proof. — Using the exact sequence (3.3 1, Lemma 3.6 (c),(d) and the fact that }i^(X,TX) = 0, we get an exact sequence 



of //-modules 

— > h' {x,tx) — >v'^ — >v — >0. 

The result follows since the characteristic polynomial of an endomorphism is multiplicative under equivariant short 
exact sequences. □ 
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Corollary 3.8. — For any matrix M of infinite order in GL(2; A), is a rigid automorphism. 



Proof. — Let M be a matrix of infi nite o rder in GL(2; A). Since H has finite index in GL(2; A), there exists an integer p 
such that M'' Hes in H. By Theorem 3.7 1 is not an eigenvalue of /^,, (and therefore of /^) . Thus /m is infinitesimally 
rigid so it is rigid. □ 



For the sake of completeness, we also provide briefly the corresponding results for cubical lattices. The same strategy 
works but the results are completely different. In the cubical case, the group G is of order 4, the map (|) (which is this 
case the multiplication by i) has 4 fixed points and 12 new other points are fixed by (|)^. After blowing up these 16 points 
and taking the quotient by G, we get a rational Kummer surface X that can also be obtained by blowing up in 1 3 
points. Besides, A is the cyclic covering of X along 10 disjoint rational curves E\,. . . ,Eio, where the multiplicity of 
the covering is 4 along £1, E2, £3, E4 and 2 along the other ten remaining Et's. We can define the subgroup K of finite 
index in SL(2; A) consisting of the elements fixing the 16 points. In the cubic case 
Therefore we get an isomorphism of /T-modules 

10 

H\x,TX)c,^n\Ei,NE,/x). 

1=1 

The first four terms correspond to the A'-module Sym^Tp.A (which is of rank 3) and the six last ones to the trivial 
/T- module C. Therefore we get 



Theorem 3.9. — For a cubic Kummer surface X the K-module (X,TX^ of infinitesimal deformations is isomorphic 

\ 4 ^ ^6 



to (Sym^y) 



o/H' {X,TX) is equal to \ Xg, 



In particular, for any element M in K, the characteristic polynomial of the endomorphism f^ 



,(x)) (. 



1 



As a corollary, we see that elements in K are not infinitesimally rigid. We can prove a slightly stronger result: 
since (0,0) is always fixed by any matrix M, H' (£1 ,N£|/x) is a GL(2; A)-submodule of H' (X,TX), it is isomorphic 
to AeiV ®Sym^V. Therefore, fu is not infinitesimally rigid as soon as (detM)^ = 1. Nevertheless it is possible to 
construct infinitesimally rigid elements in GL(2; A), but they are of lower interest than in the hexagonal case because 
they do not remain infinitesimally rigid after iteration. For instance, we have the following result: 



Theorem 3.10. — Let M be a matrix in GL{2; A) such that M 
of the endomorphism f^ o/H' (X,TX) is equal to 

4 



1 
i 



mod 2A. Then the characteristic polynomial 



(XdetM-lxSym^M-'W) {x - dttuf {x + AtiUf . 



In particular, if I is not an eigenvalue ofdetM x Sym M, then fM is infinitesimally rigid. 



4. Realization of infinitesimal deformations using divisors 

Let X be a rational surface and let / be an automorphism of X. We will present a practical method to compute the 
actions /* and /* of / on (X,TX) using divisors. This is much more delicate than for the action on the Picard group. 
In practical examples, this method is effective {see 
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4.1. 1-exceptional divisors. — We start by defining a special class of divisors which will be of crucial importance in 
our study: 

Definition 4.1. — Let X be a complex surface and let D be an effective divisor. We say that D is \-exceptional if the 
natural morphism H' {X.IX) — > H' {X,TX{D)) vanishes. 

The importance of these divisors comes from the following immediate observation: if D is 1-exceptional, then we 
have an exact sequence 

— >}^\X,TX) — >\\^{X,TX{D)) — >V]^\D,TX{D)\o) — >n\X,TX) — ^0 (4.1) 
so that if X has no nonzero holomorphic vector field, 

, }\^{DJX{D\d) 



ho(x,tx(d)) 

Let us now explain how to construct these divisors on rational surfaces. We fix two complex surfaces X and Y such 
that Y is the blowup of X at a point p. If 7i is the blowup map and E is the exceptional divisor, we have a fundamental 
exact sequence relating TX and TF via the differential of 7i which is 

O^TF ^7t*TX^l£H.ili(-l) ^0. (4.2) 

For any locally free sheaf on X, the natural pullback morphism H' (X , ^) — > H' (F, 7i* ^) is an isomorphism; this 
follows by writing down the Leray spectral sequence for the sheaf 7t* and using that R^TI* Oy =Q for j>0 {see Lemma 
|3.6| f/j where we already used this argument). Therefore, for every divisor D on X, we have a long exact sequence 

0^H''(F,TF(7i:*D)) ^H°(X,TX(D)) ^ TpX^L^ ^ (F,TF(7I*D)) ^H1(X,TX(D)) — =-0 

where 5d is the connection morphism and Lp is the fiber of Ox (D) at p. We can now prove the first basic result about 
1-exceptional divisors: 

Proposition 4.2. — IfD is a l-exceptional divisor on X, then %*D + E is l-exceptional on Y. 

Proof. — Let us start by proving that the composition of the connection morphism 5o with the natural morphism 
H' (F,TF(7t*D)) — > (F, TF(7t*D + £')) vanishes. We have a morphism of short exact sequences 

*-TF(7I*D) ^n*TX{%*D) ^Ie*^^-'^)® Oy{ti*D) ^-0 

II I 
O^TY{n*D + E) n*TX {n* D + E) -^iE*^hi~l)<» Oy{k*D + E) -^0 

where the right vertical arrow vanishes. Therefore we get a commutative diagram 

H''(i£,n^(-1) Or(7t*D)) ^° > Hi (F,TF(7t*D)) 



I 

H''(i£*ni,(-l)®0y(7i*D + £)) ^Hi(F,TF(7i*D + £)) 
and we get our claim since the left vertical arrow vanishes. We now consider the diagram 

(F,TF) ^ H' {X,TX) 



TpX^Lp^}i^{Y,TY{%*D)) ^H^{X,TX{D)) 

where the bottom line is exact. Since D is 1-exceptional, v = 0. This implies that the image of u lies in the image 
ofSo. □ 
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Let us introduce some extra notations. For any effective divisor D on a surface X, we put; 

jv{D)=n^{X,TX{D)) 
\W{D)=H^{X,TX{D)id). 
We also denote by [)(P^) the set of holomorphic vector fields on P^. 

Definition 4.3. — For any infinitely near point P in and any effective divisor D^ase on P^, we define a divisor Ds „ 



on BlsP^ as follows: if % = %n° "^^n- 



1 o 



1 7ii is the sequence of blowups defining P and E\,... ,En are the cor- 



responding exceptional divisors, then Dp = D^ where the D,'s are defined inductively by the formulas Dq — Dbase 
Ei for 1 < / < A^. If Dbase is empty, we put D^ p — Dp. 



and Di = 7t*_iA_i 



According to Proposition 



4.2 



the divisor Dp is 1 -exceptional on BlpP^. 

,2 



Note that this definition extends readily 



to a finite number of infinitely near points in 

Lemma 4.4. — For any infinitely near point P in P^ of length N, the equality dim W (Dp) — 2N holds. 

Proof. — Let X be the rational surface obtained by blowing up P. Since 7t(£>p) is a finite set in P^, V{Dp) embeds 
in [}(P^). On the other hand, it is easy to see that each holomorphic vector field on P^ yields a meromorphic vector field 
on X which is in V{Dp). Therefore V{Dp) is isomorphic to [)(P^) so that we have 



hUx.TX) 



A\mW{Dp) 

Since h^ (X,TX) - \i^iX,i:x) = 2A^- 8, we get the result. 



8-h"(X,TX 



□ 



Proposition 4.5. — Let X be a rational surface without nonzero holomorphic vector field obtained by blowing P^ 
in infinitely near points Let Dbase be an effective divisor on P^, and let y^(Z)base) be a direct factor 

o/f)(P^) in V(-Dbase)- denote by D the 1-exceptional divisor Dp^^ ^p^ . Then tliere is a natural isomorphism 
between V (Dbase) andV{D), and the associated morphism 

k 

0W(Dp.)®y^(Dbase) ^ 



;=1 



is an isomorphism. 



Proof. — For the first isomorphism, we argue exactly as in t he p roof of Lemma |4.4| For the second isomorphism, 
set K — Dp^ + ... +Dp^. Let us write down the exact sequence (4.1 
get a commutative diagram 











I)(P2) 



■V(Dba 



■W{K) 



■W{D) 



for the divisors K and D. Since V{K) = ()(P^), we 
H^{X,TX) - 



nHx.Tx) 











As y(£)base) = f)(P^) © V''^(Dbase), we obtain that W(Z)) is isomorphic to ® y■|"(Dbase)■ 
We end the section with a statement which is particularly useful in practical computations: 
Lemma 4.6. — Let D, D' be two divisors on a rational surface X such that: 

{D' is 1 — exceptional. 
The natural map from W{D) to W{D') is surjective. 

Then D is l-exceptional. 



□ 
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Proof. — We have a commutative diagram 

W{D) ^H1(X,TX) 



W{D') 



where the bottom horizontal arrow is onto. Hence the top horizontal arrow is also onto. 



□ 



4.2. Geometric bases. — In this section, we construct a basis for the vector space W{D) using coverings (the divisor D 
has been introduced in Proposition 4.5 i, we will call such a basis a geometric basis. Since TX{D)\0 is the sheaf 

TXfZ 

quotient 



a global section can be represented as a family of sections {Zi)i^j of TX(D) on open sets t/,- covering D 



such that Z, — Zj is holomorphic on Ujj for all /, j in /. 

We start with the simplest case, namely Dbase = and ^ = 1, so that D = Dp. Let us construct this basis by induction 
on the length of P. 

Let P be an infinitely near point of of length A^, let ^ be a point in one of the exceptional divisors, P' = PLi {q}, 
and put X = BlpP^ and Y = Blp, P^. Assume that we are given a covering of Dp by open sets {Ui)i(zi of X such that : 
® there exists a unique /q in / such that q G Iq; 
® the evaluation map H''(t/,„,Tt/,u) — > T^X is surjective; 

@ there exists a basis Zi, . . . ^Z2n of W{Dp), where for each a, Za is a collection of holomorphic sections {Zia)iei 
of TX{Dp) on the f/,'s such that for all / and j, the section Z,a — Zja is holomorphic on each Uia H Uja. 

Let 7t: F — > X be the blowup map and let E be the exceptional divisor Set f// = 7t^'(f/,). For any basis (vi, V2) 
of (BlpP^), thanks to @, we can choose two holomorphic vector fields Ti , T2 on Uig which extend (vi , V2). We will 
consider the vector fields Ti and T2 as sections of TY{E) on U-^. Let us now consider the sections (Zj, . . . jZj^^j) 
ofW(Dp) defined using the covering (Ul^i^j as follows : 





Z\ ... 


7' 


7' 


7' 




Zuo . . . 


■ ■ ■ Z2Nio 


Ti 


T2 


U'j j^io 


Zij ... 


■ ■ ■ Z2Nj 









Then we have the following result: 
Proposition 4.7. — The family (Z\ ^Z'^f^^^ is a basis ofW(^Dp,y 



Proof. — Thanks to Lemma 4.4 it suffices to show that this family is free. Let A-i , 
such that A-iZ'j + X2Z2 + ■ • • + X2n+2Z2n^2 — ^- Then for all j 7^ /q, 

A-iZ'iy + X,2Z2j + . . . + X2nZ2!^j 



,A,2Af+2 be complex numbers 



is holomorphic on f/J, and XiZ[-^^ + A-2Z2,| 



"2(0 " 



X2NZ2Ni^+X2N+iTi +X2N+2T2 is holomorphic on 



Therefore, X\Z\j + X2Z'2j + . . . + A,2a(Z2^^. is holomorphic on U'j also for j = /q and X2n+\P\ +X2N+2T2 (considered 
as a vector field on t/jg) vanishes a.lq. Since (Zi, . . . ,Z2n) is a basis of W {Dp) we get A-i — X2N = 0; and as (vi, V2) 
is a basis of T^^X, X2N+1 = X2N+2 =0. □ 

Remarks 4.8. — 



(i) In practical computations, f/,^ is a domain for some holomorphic coordinates (z, w), Ti ~ 3, and T2 = 3m 
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(ii) If ^' is a point on P', the covering {Ul)iei does not satisfy in general conditions ® and @. Therefore it is 
necessary to refine the covering at each blowup. 



Using Proposition 4.5 we can now easily construct a basis of W (D) . Indeed, we have already explained how to 
construct a basis of each W(^DpJ. Then it suffices to take a basis of V(£>base)^ and to pull back these meromorphic 
vector fields to the surface. 

4.3. Algebraic bases. — In practical computations, if an element of W (D) is given by local meromorphic vector fields 
on open sets of a covering, it is not a priori obvious to decompose this element in a geometric basis (as constructed 



in 5 4.2 1. In this section, we will construct another basis (we call it an algebraic basis) which solves this problem. 



We start by defining the residue morphisms. Let P = {pi,... ,pn} and let £1 , . . . ,En be the corresponding strict 
transforms of the exceptional divisors. We fix holomorphic coordinates near pi such that pi = (0,0) in these 

coordinates, and we introduce new holomorphic coordinates (m,, v,) by putting x; = m; and — UjVi. We can consider 
(m, , V,) as holomorphic coordinates on an open subset of X which contains £, except a finite number of points. 

Let (0, A,,) be a generic point of let Z be an element of W {Dpj and let Z, be a section of TX{D) which lifts Z near 
this point. Then Z, admits a Laurent expansion 

-7/ ^ (v/-A.;)"'(fl„,„3„+/7,„„3,,) . , 

Zi[Ui,Vi) = 2^2^ h holomorphic terms. 

n=lm=0 ^ 

Definition 4.9. — For any generic complex number A-,, we define the residue morphism resg^ : W[Dp) — > C" by 
the formula res£,. (Z) = {bw^bw). We also define resp : W{Dp) — > C^^ by resp = res£,.. 

Note that the definition of these residues depends on the coordinates and on the generic parameters on the exceptional 
divisors. 

Lemma 4.10. — The morphism ress is an isomorphism. 



Proof. — Let us take a geometric basis associated to the coordinates {xi,yi)\<i<N as constructed in 54.2 Then the 
matrix of resp is lower triangular by blocks, where the diagonal blocks are the 2x2 matrices whose columns are the 
vectors res£,.(8;t,) and res,E^{dy^). We have = 3„, — ^9,., and dy. = ^.^vr Therefore we obtain res£;(3;i:,.) — ( — 1, —A,) 
and res£:;(3v, ) ~ (1,0), so that the diagonal blocks are all invertible. □ 



We can now define algebraic basis in the general case. We take the notation of Proposition 4.5 In order to avoid 
cumbersome notation, we will assume for simplicity that Dbase is ureducible, and leave the general case to the reader. 

Letxi, Xr be a basis of y(Dbase)'''■ 
Le»^^Mfl 4.77. — Assume that Z^base irreducible, and let Z be an element ofW{D). Then there exist unique complex 
numbers ai (Z), . . ., (X,-(Z) such that for any generic point ^ ofD\,a^e and any lift Z ofZ near ^ , 

Z = <X\ (Z) Xi + ■ • • + ct;n(-Z) Xr + holomorphic terms 

in a neighborhood of^ . 

Proof. — The existence of such a decomposition is straightforward: the a,(Z)'s are the coefficients of Z in y(Dbase)^ 
when decomposed in a geometric basis. For the unicity, we remark that ai Xi + • • ■ + OC, Xr is holomorphic near a point 
of Z^base, it must be holomorphic on since Dbase is irreducible. Since V(£>base)^ is a direct factor of f)(P^), all the 
coefficients a, must vanish. □ 

Corollary 4.12. — With the notations of Proposition \4.5\ the map 

(resp^ , . . . ,resp^ , ai , . . . , a,) : W{D) c^N,+...+2N,+r 

is an isomorphism. 
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By definition, the associated algebraic basis of W{b) is the image of the canonical basis of c2A'i+...+2Afj.+r 
inverse isomorphism. 

5. An explicit example on blown up in 15 points 

5.1. The strategy. — Let (|): — -> P^ be the birational map given by [jc: y: z] — [xz^ : yz^'. z^]- The map (|) 
induces an isomorphism (which we still denote by (|)) between Blp^ P^ and Blp^ P^, where Pi and P2 are infinitely near 
points of P^ of length 5 centered at [1 : : 0] described in IJ] §2.1]. For any complex number a, we put 

"a 2(1 -a) 2 + 
A= -1 a+1 
1 -2 1-a 

and we consider it as an element of PGL(3; C). Then the map A(|) lifts to an automorphism \\t of the rational surface X 
obtained by blowing up the projective plane 15 times at Pi, AP2 and A(|)AP2 (see iH §3.3]). Recall that the parameter a 
is not really interesting because two different values of a correspond to linearly conjugate automorphisms, however we 
keep it in order to check some calculations in the sequel. Our aim is to compute the map acting on {X.TX). 

Let us recall the construction of Pi (resp. P2). In affine coordinates {y,z), we blow up the point (0,0)y,,; in P^ 
and we denote by E the exceptional divisor. We put y — ui, z = uivi, we blow up (0,0)„j,v[ and F is the exceptional 
divisor Then we put ui = r2S2, vi = S2, we blow up (0,0)r2,.i2 and we denote by G the exceptional divisor. Next, we 
put r2 ~ rj,Si, S2 = 53, we blow up (— I,0)r3,.s3 (resp. (1,0)^3,4-3) and H (resp. K) is the exceptional divisor. Then, we 
put r4S4 — 1, — 54 (resp. ft, — C4J4 + 1, — d^), we blow up (0,0),,,.^ (resp. (0,0)c^.^^) and we denote by L 
(resp. M) the last exceptional divisor Lastly, we put r4 — r^s^, 54 — 55 (resp. C4 — c^ds, d^ — d^Y 

Set Xi = Blp^ p2 and X2 = Blp^ P^. Let A be the strict transform of the line A = {z = 0} in X2. Then 

^(£)=£, ^(F)^K, ^(G) = G, m)^F. ^(L)=A. (5.1) 

We have Dp^^ E +2F + + + 6L wd Df^= E + 2F + AG + 5K + 6M. Let Di be the divisor on Xi given by 

Di =E + 2F + 3G + 4-H + 5L. 

By explicit calculation, the natural morphism from to W(Z)p ) is surjective so that Di is 1 -exceptional on Xi by 

Lemma |4~6| We now define a divisor D2 on X2 as follows {see Definition [43]i: 

By ( |5.1| l, we have < D2. Let J)i and 'D2 be the two 1-exceptional divisors on X given by: 
Di =Di +ADp^+A^ADp^, D2=Dp^ +AD2+A^ADp^. 

Then Si < S)2 and < 32- Therefore the morphism \|/h, acting on H' {X,TX) can be obtained as the composition 

W{Di) V. W{D2) ~ W(2)i) 



y(S)i) y(D2) y(2i) 



(5.2) 



where the inverse of last arrow is induced by the natural morphism from W (Si) to W(J)2)- To compute the strategy 
runs as follows: 

Calculations for the pair (Di,D2). 

Step 1 - Express the vectors of a geometric basis of W(Di ) in an algebraic basis. 

Step 2 - Compute (|)* : W{Di) — !> W{D2) where W{Di) is endowed with a geometric basis and W{p2) is endowed 
with an algebraic basis. 

Step 3 - Find bases of V{Di) and V{D2) whose vectors are expressed in algebraic bases of W{Di) and W{D2) 
respectively. 
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Calculations for the pair (2)1,2)2). 

Step 1 - Compute \|/* : W(2)i) — > W(2)2) where W(2)i) and W{D2) are endowed with algebraic bases. 
Step 2 - Find bases of V(2)i) and y(2)2) whose vectors are expressed in algebraic bases of W(2)i) and W(2)2) 
respectively. 

5.2. Calculations for the pair (Di ,7)2). — We give results of the calculations for the first three steps listed above. 

Step 1 - We use the coordinates {y-,z), (mi,vi), {r2,S2), (^3,53), (^4,54) and {r5,s$) to compute the algebraic and 
geometric bases on Xi . Let (f/j ) i<k5 be the covering of Di given by the following picture: 




U5 

Then a geometric basis ofW{Di) is given by the ten following vectors: 





Ui 


U2 


U3 


Ua 


U5 




dy 


dy 


dy 


dy 


dy 


62 


dz 


dz 


dz 


dz 


dz 







3m 1 


3m 1 


3m 1 


3m 1 


e4 





3vi 


3vi 


3vi 


3vi 


65 








3^2 


3r2 


3r2 


ee 








3^2 


3^2 


3*2 


ev 











3r3 


3r3 


es 











3.? 3 


3.V3 


cy 














3/4 


eio 














3*4 



- Matrix ^ - 



If X\, X2, ^3, ^4 and are generic parameters on E, F, G, H and L respectively and if (l()i<(<io is the associated 
algebraic basis, a direct calculation yields: 
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ei 


Z2 


23 


Z4 


65 


eg 


e^ 


eg 


eg 




ll 


-Xi 


1 


























[2 


-1 





























b 


2 





1 






















14 











-1 































-2X3 


1 


-X3 














k 











-2 





-1 














[7 


— X4 


K 





-3X4 





— 2X4 


1 


— X4 








ts 


—2X4 







-3 





-2 





-1 








t9 





2Xi 





-4X5 





-3X5 





-2X5 


1 


-X5 


lio 





4I5 





-4 





-3 





-2 





-1 



- Matrix - 

Step 2 - We start by fixing a basis of U°{F'^,TF'^{mA)). We divide H''(p2,Tp2(mA)) in four subspaces, the first one 
corresponding to holomorphic vector fields. 



-Vector fields of type a that span a subspace of dimension 8: 

ai=dy, a2=ydy, a-i=zdy, 04 =3^, 05 



-yd^, ae=zdz, a-]=yz^z, a%=z^dz 



-Vector fields of type b that span a subspace of dimension '"^'"^^^ ; 



l<q<m, 0<p<q + l. 



-Vector fields of type c that span a subspace of dimension "'■"^^^ : 



l<q<m, 0<p<q + l. 



-Vector fields of type d that span a subspace of dimension m: 



yP + 2 



ZP 



_1 "z 



I < p <m. 



We take for H" (P^, Tp2(mA)) ^ the subspace spanned by meromorphic vector fields of type b, c, and d (in our setting, 
m = 5). Then we use the coordinates {y,z), (mi,vi), (?"2,S2), {r3,S2,), {04, ((4) and {05, ds) on X2. If /Ji, /J2, m, H4 and 
are the generic parameters on E, F, G, K and M respectively, we consider an algebraic basis of W{D2) associated with 
these parameters consisting of ten vectors (m,)i<,<io corresponding to the exceptional divisors, and 55 vectors (25 of 
type b, 25 of type c and 5 of type d) corresponding to abasis of H'^(p2,Tp2(5A)) . 

By direct computation, we have 



<|)*3r4 = 



4j3 3y6 



1-2 



y , y 



dy- 



3/ _|_3r 



--2 

y z 
r 



y , y 



y 



-2y+^ 



-2z+' 



-2zH 
_2z2 

y 



y_ 

z 

+y 



9z, 

az, 

3z, 

^ az, 



^^dz = 2 
<l)*3vi = 

^*^S2 = 



y__y_ 
^ z^ 

.5 



av+ 1 + 



2y^ 



2^8,,+ U + 2 



z 



'y + '^y^z, 
..2- 



Uy 



2^ Hz- 
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Therefore, the matrix of (|), : W{Di) — > W{D2) is: 





<|)*ei 








(l)*e5 


*l)*e6 


c|)*e7 


c|)*e8 


(l)*e9 


(l)*eio 


mi 




1 


























xn2 


-1 




























m3 


2 





1 


-m 


-2/J2 







m 








nT4 











-1 


-2 







1 








ms 











2^3(^3-1) 




/'3(W-1) 














me 











4/i3-2 


2(/i3-l) 


2/i3-l 














XtlJ 






A'4 


A'4 




















m 




-H 


2^4 


1 




















my 


































































-4^3.2 


2/j4,3 


-2^3.2 




-2^2,1 


^2,1 




2/72,1 






H0(p2,Tp2(5A))^ 


-3C2,1 

3d4 


C3,2 
-2d5 


-2C2,1 

d4 


2d3 


d3 


^^3 


di 




di 





- Matrix ^ - 



Step 3 - For typographical reasons, we will write the vectors of y(Di) and V{D2) as lines (and not as columns). 





u 


h 


13 


k 


(5 


k 


h 


k 


[9 


lio 


ai 




-1 


2 











— A-4 


—2X4 








ai 




















-3 











0-3 


























-3 





a4 


1 

















K 


3Xi 


2Xi 


4?15 


as 








-^2 


-1 


-2A.3 


-2 


-3X4 


-3 


-4X5 


-4 






















2 











ai 
































as 










































- Matrix 












mi 


m2 


m3 


m4 


ms 


me 


m^ 


mg 


m9 


triio 


a\ 




-1 


2 













-2^X4 








02 




















3 











03 


























3 





an 


1 



















-3^*4 




4//5 


as 








-m 


-1 


-2/J3 


-2 


-3/J4 


-3 


-4a/5 


-4 






















-2 











































as 

































- Matrix - 
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mi 


m2 


m3 


m4 




mg 


rriT 


trig 




TTllO 






















2a'4 


• 4 


3a/^ 


9A'i 


bi,i 


-1 

















A*! 


r q. 


~j 


-4p5 


bi,! 








2/J2 


2 


3/J3 


3 


4/J4 


4 


Sa's 


5 


bo,2 






























^1,2 




















-3a'4 

' 4 


' 4 








bi.i 


-at' 


at' 














—/A 


-4a'1 








bi.2 
































bQ3 




















52^/^ 


572^1" 






bi,3 




















r 4 


135^^ 




60Ai5 


^2,3 




















4a'1 


28/i$ 

("4 




-20/iS 


^3.3 




"1 


















2ivk 


6ui 


^4,3 
































boA 






















-4760^/^^ 








b\A 




















-9lA'i^ 










b2,A 






















-220^/^ 








b3,4 






















-40^/^ 








b4,4 




Sat' 














~4 


-fiA/l 








bsA 
































bo,5 




















2394^1' 


40698A/f 


429^1 




^1,5 




















612^f 






-1320^^ 


b2.,5 




















140A/f 


1820^/^^ 




420Ai^ 


bi,5 




















30A/V 




-21^1 




bAS 




















6a'4 




7A'^ 


35a/^ 


bs.s 


-AT^ 
















/4 




-2^t 




b6,5 

































- Matrix - 
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mi 


m2 


m3 


m4 


ms 


me 




mg 


mg 




co,i 




















4^*4 


24^^ 








Cl,l 


at' 

' 1 


-at' 














r 4 










C2,l 
































C0,2 


























-904 


Cl,2 




















r 4 


-35^5 

r 4 




244 


C2 2 




1 
















— 5wl 

r" 4 


-2u\ 




C3,2 
































C03 






















2184py 








C1.3 






















330^^ 








C2,3 




















(vl 










C3,3 


at' 


" 1 
















r*4 








C4,3 
































C0,4 




















-1428^/^^ 




330^^ 


2640a/^ 


Cl,4 






















-3185^/^^ 




-6724 


C2,4 






















-462^1" 




168^^ 


C3,4 






















-63/i« 


-s4 


-404 


C4,4 


AT^ 


-4aT' 














-A 




24 


s4 


C5,4 
































C0,5 




















11704A/f 


210672/7^ 








Cl,5 




















1938A/f 


31008/if 








C2,5 




















320^/^* 


4480A/f 








C3,5 




















52^/^ 


624A/V 








C4,5 






















80^/^ 








C5,5 


A^r^ 
















A'4 










C6,5 

































- Matrix '.yV^ - 





mi 


m2 


m3 


m4 


11^5 


mg 




mg 


mg 


mio 


d\ 


























1 





ih 










2l.i2 








1 











ds 
















2^3 


2(04 


2 


2a'5 


2 


C?4 




















^'i 


2/J4 








d5 


























A*^ 


2/i5 



- Matrix 'j^i - 



Although we won't need it, it is easy to verify that 

(])*(ai) =ai -4fc3,2-3c2,i +3<i4, (|)*(a2) =(22-3^2, ^*{a^) = a^-3di, 

(|)H.(fl4) =fl4- 2^4,3 +C3,2-2<i5, (])* (fls ) = as + 2<i3 , (|),(a6) =a6 + 2(i2, (5.3) 

(|)H.(fl7) = fly, (|)H.(a8) = fl8- 

5.3. Calculations for the pair (2)1,2)2). — In this part we provide the last two steps of the calculations. 

Step 1 - We have isomorphisms 

W{T)i) ~W (Di) ®W {A ■ DjsJ ®W {AfA ■ Dp^) , W(D2) ~ W(Dp^) ® W(AZ)2) ® W(A/A -DpJ . 

We transport the bases of W {D2 ) and W {Dp^ ) by A and A^A. Therefore, if we take algebraic bases of W (Dp^ ) , W (Dp^ ) 
and W{D2), the matrix of \|/h. : W{Di) — > W{D2) has the form 
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A(|)-(l,)i<i<io 


A(|)- (Am,)i<i<io 


A(|)- (A(|)Am,)i<,<io 


ili)l<i<W 


OlOxlO 


OlOxlO 




(Am,)i<Kio 
H''(P^TP2(5AA))''' 




055x10 


055x10 


(A(|)Am,)i<(<io 


OlOxlO 


Idioxio 


OlOxlO 



The matrix ^ is the 10 x 10 matrix given exphcitly by 



-4:^1 - 1 


qi.2 


























-4 


4/ii-l 


























8-25:^2 


/ii(25?02-8) 


6^2-1 


93,4 




















-25 


25/Ji 


6 


1 -6(U2 




















-50^3 


8 + 50^3/ii 


4(l+3?^3) 


-4/i2(l + 3X3) 


-1 


i"3+^3 














-50 


50/ii 


12 


-12/i2 





1 














17,1 


17,2 


18:^4-42 


/i2 (42 -18^4) 


-10 


10;U3 - 10 


-1 


;U4 + ^4 












18 


-18/i2 











1 








19,1 


19,2 


19,3 


99,4 


-110 


110Ai3-62 


-12 


12/M 


-1 


99,10 


^10,1 


1W,2 


24 


-24/i2 

















1 



where 

91,2 =<"!(! +4Xi)-Xi, <73,4=A'2(l-6?i2)+?i2, <?7,i = 2a + 200 + ?14(3a/^ - - 75 - 4X4) + 2a/4, 
97,2 = 40 + A/1 (75A.4 - 2a - 200 + 4X1 - 2/4 ~ 3^/^X4 + X|) - {lx^ + , q%.i = 3(^/1 - X^) - 75 - 8X4, 
98,2 = li\ (75 + 8X4 + 3(X| -A'4)) - 2(^(4 +X4), <79,i = 672 -4^/5 (X5 + 1) +24(a/^ + a) - 2(X^ -\ix\) - IOOX5, 
99,2 = A'l (4^5(1 + >^5) - 24(a/1 + a) + 2(^^ + X^) + IOOX5 - 672) + 320 - 2a - 12^(4, ^9,3 = 24X5 - 92 - 4a, 
99,4=^*2 (92 + 4a -24X5), ^9,10 =/^5 + 1+^5, 910,1 = -4(25+ X5+Ai5), ^10,2 = 4a/i(25 +X5 +//5). 

5rep2-Wehavey(2)i) = ()(p2) andy(2)2) =H0(p2,Tp2(5AA)). Let us decompose a basis ofy(S)i) in an algebraic 
basis of W(2)i): 





(«()l<K8 


(l/)l<i<10 




(Am,)i<Kio 




(A(l)Ami)i<Kio 





- Mam'x >i - 



To compute quickly the matrices and ^a, we make the following remark (which can be proved by an easy 
computation): for any holomorphic vector field defined in a neighborhood (0,0), if we lift it as a section of W(Dp^) 
or ly(Dp^), this section depends only on the seven Taylor components 3y, y^y, zdy ,3^, yd^, z9z, y^^z of the vector field 
at (0,0). The corresponding sections in algebraic bases are given by the tables 
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dy 








yd. 


zd. 


y'd. 


ll 










1 











il 


-1 




















h 


2 



















u 














-1 








h 
















-2A.3 








k 














-2 








h 


— A,4 


-3 









-3A.4 


2 





k 


— 2A.4 










-3 








h 








-3 


2X1 




-4X5 





-2 


ho 















-4 


















- Matrix 3fi 
















dy 

} 






5z 




yd. 




zd. 




y'd. 


m 










1 











m.2 


-1 




















m3 


2 











-A'2 








m4 














-1 






















-2/J3 








mg 














-2 








mv 


-A'4 


3 





-A'4 


-3/J4 


-2 





mg 


-2/J4 










-3 








mg 








3 


24 


-4,^5 





-2 













4/i5 


-4 









- Matrix 2^2 - 



Therefore, in order to compute J^'^ and J^"^, it suffices to extract the seven aforementioned Taylor components of 
the vector field A~^a, (resp. (A(|)A)~^a,) at (0,0). Then we multiply the resulting vector by the matrix ^1 (resp. ^) 
and we obtain (resp. J^"^ ■ We won't give the exact expressions of and because of lack of space. 

We now deal with Vi^q). We decompose a basis of vectors of Vi^i) in an algebraic basis of W(S2). We get the 
following matrix 





(Aa,)i<,<8 






{Adp)i<p<m 


(t,)l<K10 










(Am,)i<,<io 

(^^p,9)o<p<9+l<m+l 
(^Cp_^)o</)<9+l<m+l 
[Adp)\<p<m 


^a 
025 x8 
025 x8 

05x8 


W25x25 
025x25 
05x25 


025x25 

Id25x25 
05x25 


025 x 25 
025 x 25 
Id5x5 


(A(|)Am,)i<Kio 




^' 







- Matrix >2 - 



The matrices (resp. =/^a)' (resp. ^b), (resp. ./K") and (resp. ^^') appearing in >2 are 
computed in the same way as (resp. ^!^) . 

5.4. The result. — We can now state and prove the following result. 
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Theorem 5.1. — Let X be the rational surface obtained by blowing up the projective plane 15 times at the infinitely 
near points P\, AP2 andA^AP2 and l|/ be the lift ofA^ as an automorphism ofX. Then the characteristic polynomial 
of the map \|/* acting on the space H' (X,TX^ of infinitesimal deformations ofX is 

Q^{x) = {x^ + 3x+l) (x^ + lSx+l) (x^-Vx+l) (x^+x+l) {x - if {x + l)\x^ - x + l)"^ . 
Besides, there is only one nontrivial Jordan block, which is a 2 x 2 Jordan block attached with the eigenvalue —1. 
Proof. — Let tfi be the subspace of W(S)i) of dimension 22 defined by 

{h, hi I3, 14, [5, ki h, kj k, ho, 1 
Ami, Am2, Ams, Am4, Ams, Ame, Amv, / ■ 
A(|)Ami, A(|)Am2, A(|)Am3, A(|)Am4, A(|)Am5 J 

A direct calculation shows that S'l is a direct factor of V(£'i) in W{'I)i). If (S2 is the image of (fi in W{'I>2) by the 
natural injection, then £2 is a direct factor of V (2I2) in W{'I)2). 

Therefore, the composition of the morphisms ( |5.2| i can be expressed as 

A ® ~ W(£>2) ~ ^2 © V(S)2) — ^ <?2^<^i ■ 

The matrix of /' (resp. p) can be computed using Yi (resp. Y2) and the matrix of \(/* is ^V. □ 

Corollary 5.2. — If (X, f) is a holomorphic family of rational surfaces endowed with an automorphism that lifts the 
pair {X,\\t), then m(X) < 2. 

It is an interesting problem to know if this bound for m(X) is sharp. 
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> restart: 

> with (Linear Algebra) : 



6.1 First definitions 

> restart: 

> with(LinearAlgebra) : 

The function f is the function of the text given in the affine coordinates (y, z). 

> fl:=(y,z)->(y*z'2)/(z-2+y-3): 

> f2:=(y,z)->(z-3)/(z-2+y-3) : 

> f :=(y,z)->(fl(y,z) ,f2(y,z)) : 

The inverse of f is denoted by g. 

> gl:=(y,z)->y*z-2/(z~2-y-3) : 

> g2:=(y,z)->z"3/(z"2-y'3) : 

> g: = (y,z)->(gl(y,z) ,g2(y,z)) : 

We define the action of the matrix A in the coordinates (y, z), its inverse is B. 

> A:=(y,z)->((-l+(l+alpha)*z)/(alpha+2*(l-alpha)*y+(2+alpha-alpha~2)*z) , (l-2*y+(l-alpha) * 
z) / (alpha+2* ( 1-alpha) *y+ (2+alpha-alpha"2) *z) ) : 

> B:=(y,z)->((l-y-(l+alpha)*z)/(l+alpha+(alpha-3)*y+(l-alpha~2)*z) , (l+y+( 1-alpha) *z)/(l+alpha+ 
(alpha-3)*y+(l-alpha"2)*z)) : 

We define m = Af Af A and n = BgBgB. 

> ml : =(y, z)->op(l, [A(f(A(f(A(y,z) ))))]) : 

> m2:=(y,z)->op(2, [A(f (A(f(A(y,z) ))))]) : 

> m:=(y,z)->(ml(y,z) , m2(y,z)): 

> nl: = (y,z)->op(l, [B(g(B(g(B(y,z) ))))]) : 

> n2: = (y,z)->op(2, [B(g(B(g(B(y,z)))))] ) : 

> n:=(y,z)->(nl(y,z) , n2(y,z)): 



6.2 The procedures to compute direct images of vector fields 
6.2.1 The procedures chcimp 

The procedure champ 1 associates to a vector field Z the coordinates of the Taylor expansion of the vector field 
A-t Z in the following order: dy, ydy, zdy, dz, ydz, zdz, dz- The variable of the procedure is not the vector 
field itself but a curve t i— >■ Mt whose the derivate at is Z. For example, if Z = u(y, z)dy + v{y, z)dz, we can 
take Mt{y,z) = {u{y, z) + ty,v{y, z) +tz). 

> chcimpl : =proc (M) 

> local prod, res, resl,res2,u,v: 

> global resul ,resu2,w: 

> prod:=unapply([(A@(M@B)) (y,z)] ,t) : 

> res : =simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply(op(l ,res) , (y ,z) ) : 

> res2:=unapply(op(2,res) , (y,z)) : 

> resul : = [subs ({y=0 ,z=0} , resl (y , z) ) , subs ({y=0 , z=0} , dif f (resl (y , z) , y) ) , subs ({y=0 , z=0} , dif f ( 
resl (y,z) ,z) )] : 

> resu2 : = [subs ({y=0 , z=0} , res2 (y , z) ) , subs ({y=0 , z=0} , dif f (res2 (y , z) , y) ) , subs ({y=0 , z=0} , dif f ( 
res2(y,z) ,z)) ,subs({y=0,z=0},diff (diff (res2(y,z) ,y) ,y)/2)] : 

> u:=op(resul) :v:=op(resu2) : 

> w : = [u , v] : 

> RETURN ( convert (w, Vector)) 

> end: 

The procedure champ2 associates to a vector field Z the coordinates of the Taylor expansion of the vector field 
{(pA)^^ Z in the following order: dy, ydy, zdy, dz, ydz, zdz, y^ dz- 

> chainp2 : =proc (M) 

> local prod, res, resl, res2,u,v: 

> global resul ,resu2,w: 

> prod:=unapply([(B(S(g@(M@(f(aA)))) (y,z)] ,t) : 

> res :=simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply (op(l ,res) , (y , z) ) : 

> res2:=unapply(op(2,res) , (y,z)) : 

> resul : = [subs ({y=0 ,z=0} , resl (y , z) ) , subs ({y=0 , z=0} , diff (resl (y , z) , y) ) , subs ({y=0 , z=0} , diff ( 
resl (y,z) ,z) )] : 

> resu2 : = [subs ({y=0 ,z=0} ,res2 (y , z) ) , subs ({y=0 , z=0} , diff (res2 (y , z) , y) ) , subs ({y=0 , z=0} , diff ( 
res2(y,z) ,z)) , subs({y=0,z=0}, diff (diff (res2(y,z) ,y) ,y)/2)] : 

> u:=op(resul) :v:=op(resu2) : 

> w : = [u , v] : 

> RETURN ( convert (w, Vector)) 

> end: 

The procedure champS associates to a vector field Z the coordinates of the Taylor expansion of the vector field 
{A(j)A)-^ (Z) in the following order: dy, ydy, zdy, dz, ydz, zdz, y"^ dz- 

> champS : =proc (M) 

> local prod, res, resl, res2: 

> global resul, resu2: 

> prod:=unapply([(B(a(g@(B(3(M@(A@(f@A)))))) (y,z)] ,t) : 

> res : =simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply(op(l ,res) , (y ,z) ) : 

> res2:=unapply(op(2,res) , (y,z)) : 

> resul : = [subs ({y=0 ,z=0} , resl (y , z) ) , subs ({y=0 , z=0} , diff (resl (y , z) , y) ) , subs ({y=0 , z=0} , diff ( 
resl(y,z) ,z))] : 

> resu2 : = [subs ({y=0 , z=0} ,res2 (y , z) ) , subs ({y=0 , z=0} , diff (res2 (y , z) , y) ) , subs ({y=0 , z=0} , diff ( 
res2(y,z) ,z)) , subs({y=0,z=0}, diff (diff (res2(y,z) ,y) ,y)/2)] : 

> RETURN (Trainspose (convert ( [seq(factor(op(k, resul)) ,k=l. . nops (resul) ) , seq(f actor (op(k,resu2) 
) ,k=l. .nops(resu2))] ,Matrix))) : 

> end: 

The procedure chainp4 associates to a vector field Z the coordinates of the Taylor expansion of the vector field 
A^^ (Z) in the following order: dy, ydy, zdy, dz, ydz, zdz, y^dz- 

> chainp4:=proc(M) 

> local prod, res, resl, res2,u,v: 

> global resul ,resu2,w: 

> prod:=unapply([(B@(M@A)) (y,z)] ,t) : 

> res : =simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply(op(l ,res) , (y ,z) ) : 

> res2:=unapply(op(2,res) , (y,z)) : 



> resul : = [subs ({y=0 ,z=0} , res 1 (y ,z) ) , subs({y=0 ,z=0} , dif f (resl (y ,z) ,y) ) , subs ({y=0 , z=0} ,dif f ( 
resKy.z) ,z))] : 

> resu2 : = [subs ({y=0 , z=0} , res2 (y , z) ) , subs ({y=0 , z=0} , dif f (res2 (y , z) , y) ) , subs ({y=0 , z=0} , dif f ( 
res2(y,z) ,z)) , subs ({y=0 ,z=0} , dif f (dif f (res2(y,z) ,y) ,y)/2)] : 

> u:=op(resul) :v:=op(resu2) : 

> w : = [u , v] : 

> RETURN ( convert (w, Vector)) 

> end: 

The procedure chaimpS associates to a vector field Z the coordinates of the Taylor expansion of the vector field 
{A(^A(^A)~^ (Z) in the following order: dy, ydy, zdy, d^, yd^, zd^, d^. 

> champS : =proc (M) 

> local prod, res, resl, res2,u,v: 

> global resul ,resu2,w: 

> prod:=unapply([(m(S(M@n)) (y,z)] ,t) : 

> res:=(subs(t=0,diff (prod(t) ,t))) : 

> resl : =unapply(op(l ,res) , (y ,z) ) : 

> res2:=unapply(op(2,res) , (y,z)) : 

> resul : = [subs ({y=0 ,z=0} , resl (y ,z) ) , subs ({y=0 ,z=0} , dif f (resl (y ,z) ,y) ) , subs ({y=0 , z=0} ,dif f ( 
resl (y,z) ,z) )] : 

> resu2 : = [subs ({y=0 ,z=0} ,res2 (y , z) ) , subs ({y=0 , z=0} , dif f (res2 (y , z) , y) ) , subs ({y=0 , z=0} , dif f ( 
res2(y,z) ,z)) ,subs({y=0,z=0},diff (diff (res2(y,z) ,y) ,y)/2)] : 

> u:=op(resul) :v:=op(resu2) : 

> w : = [u , v] : 

> RETURN ( convert (w, Vector)) 

> end: 

This procedure is too heavy for Maple, we will not use it. To replace it, we will do successive Taylor expansions 
at each step of blow-ups. We deal with it in the next section. 



6.2.2 The procedures devlim and extract 

We compute the points where we take the Taylor expansions. 

> A(0,0); 

> simplify(op(l, [(AOfOA) (0,0)])) , simplify (op(2 , [(AOfOA) (0,0)])) ; 

> simplify(op(l, [(ASf OAOf OA) (0,0)])) , simplify (op (2, [(AOfOAOfOA) (0,0)] )) ; 

0, 

The procedure devliml associates to a vector field Z in coordinates (y, z) (still given by the derivative at of a 
curve Mt) the Taylor expansion of A^Z of order 2 at a point {p, q). The result is viewed as a polynomial vector 
field in a neighborhood of (p, q) and the procedure returns a curve Nt associated to this vector field. 

> devliml : =proc (M,p , q) 

> local prod, res , resl ,res2 , resul ,resu2 ,N,MMM: 

> MMM:=(y,z)->(op(M(y,z))) : 

> prod:=unapply([(A(3(MMM@B)) (y,z)] ,t) : 

> res : =simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply(op(l ,res) , (y ,z) ) : 

> res2:=unapply(op(2,res) , (y,z)) : 

> resul : =simplif y (mtaylor (resl (y ,z) , [y=p,z=q] , 3) ) : 

> resu2:=simplify(mtaylor(res2(y,z) , [y=p,z=q] ,3)) : 

> N: =unapply ( [y+t*resul ,z+t*resu2] , (y , z) ) : 

> RETURN (N): 

> end: 

The procedure devlim2 is similar to devliml but computes the Taylor expansion of (A(f>)^Z instead of the 
Taylor expansion of A^,Z. 

> devlim2:=proc(M,p,q) 

> local prod, res , resl ,res2 , resul ,resu2 ,N,MMM: 

> MMM: = (y,z)->(op(M(y,z))) : 

> prod:=unapply([(A@f@(MMM@g@B)) (y,z)] ,t) : 

> res : =simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply (op(l ,res) , (y , z) ) : 

> res2:=unapply(op(2,res) , (y,z)) : 

> resul : =simplify (mtaylor (resl (y,z) , [y=p,z=q] , 3) ) : 



> resu2:=simplify(iiitaylor(res2(y,z) , [y=p,z=q] ,3)) : 

> N: =unapply ( [y+t*resul ,z+t*resu2] , (y , z) ) : 

> RETURN (N): 

> end: 

The procedure extract associates to a vector field given by a curve Nt{y,z) the components in dy, ydy, zdy, 
dz, ydz, zdz and y^dz in the Taylor expansion of the associated vector field at (0,0). 

> extract : =proc (M) 

> local prod, res, resl,res2,w: 

> prod:=unapply(M(y,z) ,t) : 

> res:=simplify(subs(t=0,diff (prod(t) ,t))) : 

> resl : =unapply(op(l ,res) , (y ,z) ) : 

> res2:=unapply(op(2,res) , (y,z)) : 

> w[l] :=simplify(resl(0,0)) : 

> w[2] : =simplify (subs ({y=0,z=0},diff (resl (y,z) ,y))) : 

> w [3] : =simplif y(subs({y=0 ,z=0} , dif f (resl (y ,z) , z) ) ) : 

> w[4] :=simplify(res2(0,0)) : 

> w[5] :=simplify(subs({y=0,z=0},diff (res2(y,z) ,y))) : 

> w [6] : =simplif y(subs({y=0 ,z=0} ,dif f (res2(y ,z) ,z) ) ) : 

> w[7] :=simplify(subs({y=0,z=0},diff (res2(y,z) ,y$2)/2)) : 

> w: = [w[l] ,w[2] ,w[3] ,w[4] ,w[5] ,w[6] ,w[7]] : 

> RETURN (convert (w, Vector) ) : 

> end: 

The function psij^ is the projection of the blow-up of {0,0)y^z, and etai is its inverse. 

> psil: = (u,v)->(u,u*v) : 

> etal: = (y,z)->(y,z/y) : 

The function psig is the projection of the blow-up of (— 1/a, l/a)y,z and eta2 is its inverse. 

> psi2:=(u,v)->(u-l/alpha,u*v+l/alpha) : 

> eta2: = (y,z)->(y+l/alpha, (z-l/alpha) / (y+l/alpha) ) : 

The function psig is the projection of the blow-up of (l/a, l/a)y^z and etaa is its inverse. 

> psi3:=(u,v)->(u+l/alpha, u*v+l/alpha) : 

> eta3:=(y,z)->(y-l/alpha, (z-l/alpha) / (y-l/alpha) ) : 

The procedure uldevliml associates to a vector field Z on blown up at (0,0)j,.^ written in the coordinates 
(ui, vi) the Taylor expansion of A^Z of order 3 at a point (p, q). The result is considered as a polynomial vector 
field in a neighborhood of (p, q) in the coordinates (wi, vi) which are this time the coordinates on blown up 
at {—-^,-^)y^, and the procedure gives a curve Nt associated to this vector field. 

> uldevliml :=proc(M,p,q) 

> local prod, res , resl ,res2 ,resul ,resu2 ,N,MMM: 

> MMM:=(u,v)->(op(M(u,v))) : 

> prod:=unapply([(eta2@(A(a(psil@MMM(aetal@(B@psi2)))) (u,v)] ,t) : 

> res : =simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply (simplify (op (1 ,res) ) , (u, v) ) : 

> res2:=unapply(simplify(op(2,res)) , (u,v)) : 

> resul : =simplif y (mtaylor (resl (u, v) , [u=p, v=q] ,4) ) : 

> resu2:=simplify (mtaylor (res2(u,v) , [u=p,v=q] ,4)) : 

> N : =unapply ( [u+t*resul , v+t*resu2] , (u, v) ) : 

> RETURN (N): 

> end: 

We compute the points where we take the Taylor expansions. 

> limit (op (1 , [(eta2(3(A@psil)) (u,v)] ) ,{u=0,v=0}) , limit (op (2, [(eta2@(A@psil) ) (u,v)] ) ,{u=0, 
v=0}); 

0, -{1-a)-' 

> limit (op (1 , [(eta3(a(A<af @A)(apsil) (u,v)] ) ,{u=0,v=0}) , limit (op (2, [(eta3<a(A@f (aA)(apsil) (u,v) 
]),{u=0,v=0}); 

0, -(-l-a)-i 

> limit(op(l , [(etal(a(A@f QAOf (aA)(apsil) (u,v)] ) ,{u=0,v=0}) , limit(op(2, [(etal(a(A(af QAQf (aA)Qpsil) 
(u,v)]) ,{u=0,v=0}) ; 

0, 



The procedure uldevlim2 associates to a vector field Z on the blow up of at expressed in 

the coordinates (u\,v\) the Taylor expansion of Z of order 3 at a point (p, g). The result is viewed 

as a polynomial vector field in a neighborhood of (p, in the coordinates {u\,v\) which are this time some 
coordinates on blown up at ^, and the procedure returns a curve Ni associated to this vector field. 

> uldevlim2:=proc(M,p,q) 

> local prod, res ,resl ,res2 ,resul ,resu2 ,N,MMM: 

> MMM:=(u,v)->(op(M(u,v))) : 

> prod:=unapply([(eta3@(A(Sf(S(psi2@MMM(3eta2@g@B@psi3))) (u,v)] ,t) : 

> res : =simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply (simplify (op (1 ,res) ) , (u, v) ) : 

> res2:=unapply(simplify(op(2,res)) , (u,v)) : 

> resul : =simplif y (mtaylor (resKu, v) , [u=p, v=q] ,4) ) : 

> resu2:=simplify(mtaylor(res2(u,v) , [u=p,v=q] ,4)) : 

> N: =unapply ( [u+t*resul , v+t*resu2] , (u, v) ) : 

> RETURN (N): 

> end: 

The procedure uldevlimS associates to a vector field Z on the blow up of P^ at ]^y^ expressed in the 
coordinates (u\^vi) the Taylor expansion of (Ai/))* Z of order 3 at a point {p,q). The result is viewed as a 
polynomial vector field in a neighborhood of {p,q) in the coordinates {ui,vi) which are this time coordinates 
on P^ blown up at {0,0)y^z, and the procedure gives a curve Nt associated to this vector field. 

> uldevlimS :=proc(M,p,q) 

> local prod, res , resl ,res2 , resul ,resu2 ,N,MMM: 

> MMM:=(u,v)->(op(M(u,v))) : 

> prod:=unapply([(etal@(A@f@(psi3@MMM@eta3@g@B(apsil))) (u,v)] ,t) : 

> res : =simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply (simplify (op (1 ,res) ) , (u, v) ) : 

> res2:=unapply(simplify(op(2,res)) , (u,v)) : 

> resul : =simplify (mtaylor (resl (u,v) , [u=p, v=q] ,4) ) : 

> resu2:=simplify (mtaylor (res2(u,v) , [u=p,v=q] ,4)) : 

> N : =unapply ( [u+t*resul , v+t*resu2] , (u, v) ) : 

> RETURN (N): 

> end: 

The procedure ulextract associates to a vector field given by a curve Nt{ui, vi) the components in d^, ui d^, 
Vidui, uiVidui, v^du^, vfdu^, i9„j, Vid^-^ and vfdy-^ of the Taylor expansion of the associated vector 

field at (0,0). 

> ulextract : =proc (M) 

> local prod, res, resl, res2,w: 

> prod: =unapply(M(u, v) ,t) : 

> res : =simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply(op(l ,res) , (u, v) ) : 

> res2:=unapply(op(2,res) , (u,v)) : 

> w[l] :=simplify(resl(0,0)) : 

> w[2] :=simplify(subs({u=0,v=0},diff (resl(u,v) ,u))) : 

> w [3] : =simplif y(subs({u=0 , v=0} ,dif f (resl (u, v) , v) ) ) : 

> w [4] : =simplif y(subs({u=0 , v=0} ,dif f (resl (u, v) ,u, v) ) ) : 

> w[5] :=simplify(subs({u=0,v=0},diff (resl(u,v) ,v$2))/2) : 

> w[6] :=simplify(subs({u=0,v=0},diff (resl(u,v) ,v$3))/3) : 

> w[7] :=simplify(res2(0,0)) : 

> w[8] :=simplify(subs({u=0,v=0},diff (res2(u,v) ,u))) : 

> w[9] :=simplify(subs({u=0,v=0},diff (res2(u,v) ,v))) : 

> w[10] :=simplify(subs({u=0,v=0},diff (res2(u,v) ,v$2)/2)) : 

> w: = [w[l] ,w[2] ,w[3] ,w[4] ,w[5] ,w[6] ,w[7] ,w[8] ,w[9] ,w[10]] : 

> RETURN(convert(w,Vector)) : 

> end: 

The function ppsi-^ is the projection of the blow up of (0;0)j^_2 and (0,0)„j_^j, and eetai is its inverse. 

> ppsil : =(r , s)->(r*s ,r*s'2) : 

> eetal : =(y ,z)->(y"2/z ,z/y) : 

The function ppsij is the projection of the blow up of (— ^, ^) and (0, ) , and eeta2 is its inverse. 

> ppsi2:=(r,s)->(r*s-l/alpha, r*s*(s+l/(alpha-l))+l/alpha) : 

> eeta2:=(y,z)->((y+l/alpha)/((z-l/alpha)/(y+l/alpha)-l/(alpha-l)) , (z-l/alpha)/(y+l/alpha) 
-l/(alpha-l)) : 



The function ppsig is the projection of the blow up of ( — , — ) and (0, ^|-r ) , and eetaa is its inverse. 

> ppsi3:=(r,s)->(r*s+l/alpha, r*s*(s+l/(alpha+l))+l/alpha) : 

> eeta3:=(y,z)->((y-l/alpha)/((z-l/alpha)/(y-l/alpha)-l/(alpha+l)) , (z-l/alpha)/(y-l/alpha) 
-l/(alpha+l)) : 

We compute the points where we take the Taylor expansions. 

> limit (op (1 , [(eeta2(§(A@ppsil)) (r ,s)] ) ,{r=0,s=0}) , limit (op (2, [(eeta2@(A@ppsil)) (r ,s)] ) ,{r=0, 
s=0}) ; 

0, 

> limit(op(l, [(eeta3a(A(afaA)(appsil) (r,s)] ) ,{r=0,s=0}) , limit(op(2, [(eeta3@(A@f @A)<appsil) ( 
r,s)]),{r=0,s=0}); 

0, 

> limit(op(l, [(eetal(a(A(af(aA(af (aA)@ppsil) (r,s)]) ,{r=0,s=0}) , limit(op(2, [(eetal(a(A@f OAOf OA) 
Qppsil) (r,s)]) ,{r=0,s=0}) ; 

0, 

The procedure r2devliml associates to a vector field Z on the blow up of at (0,0)y.z and (0, 0)ui,vi expressed 
in coordinates (r2,S2) the Taylor expansion of A^Z of order 2 at a point {p,q). The result is viewed as a 
polynomial vector field in a neighborhood of (p, q) in the coordinates (r2, S2) which are this time coordinates on 

blown up at (— ^, ^) and (0, ) , and the procedure returns a curve Nt associated to this vector 

field. 

> r2devliml:=proc(M,p,q) 

> local prod, res ,resl ,res2 ,resul ,resu2 ,N,MMM: 

> MMM:=(r,s)->(op(M(r,s))) : 

> prod:=unapply([(eeta2@(A(a(ppsil@MMM(aeetal@(B(appsi2)))) (r,s)] ,t) : 

> res : =simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply (simplify (op (1 ,res) ) , (r , s) ) : 

> res2:=unapply(simplify(op(2,res)) , (r,s)) : 

> resul : =simplif y (mtaylor (resKr , s) , [r=p, s=q] , 3) ) : 

> resu2 : =simplify (mtaylor (res2(r , s) , [r=p,s=q] ,3)) : 

> N : =unapply ( [r+t*resul , s+t*resu2] , (r , s) ) : 

> RETURN (N): 

> end: 

The procedure r2devlim2 associates to a vector field Z on the blow up of P'^ at (— ^, ^) and ( 0, ) 

expressed in coordinates (7'2,S2) the Taylor expansion of Z of order 2 at a point {p,q). The result is 

viewed as a polynomial vector field in a neighborhood of (p, q) in the coordinates (r2, S2) which are this time 

the coordinates of P^ blown up at ^) and (O, jjxx ) , and the procedure returns a curve Nt associated 

to this vector field. 

> r2devlim2:=proc(M,p,q) 

> local prod, res, resl, res2, resul, resu2,N,MMM: 

> MMM:=(r,s)->(op(M(r,s))) : 

> prod:=unapply([(eeta3@(A(3f(a(ppsi2(aMMM@eeta2@g@B(3ppsi3))) (r,s)] ,t) : 

> res : =simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply (simplify (op (1 ,res) ) , (r , s) ) : 

> res2:=unapply(simplify(op(2,res)) , (r,s)) : 

> resul : =simplify (mtaylor (resl (r , s) , [r=p, s=q] , 3) ) : 

> resu2 : =simplify (mtaylor (res2(r , s) , [r=p,s=q] ,3)) : 

> N : =unapply ( [r+t*resul , s+t*resu2] , (r , s) ) : 

> RETURN (N): 

> end: 

The procedure r2devlim3 associates to a vector field Z on the blow up of P^ at ^) and (0, ) 

expressed in coordinates (r2,S2) the Taylor expansion of {A(f))^,Z of order 2 at a point {p,q). The result is 
viewed as a polynomial vector field in a neighborhood of {p,q) in the coordinates (r2,S2) which are this time 
the coordinates of P^ blown up at (0,0)^,^ and {0,0)ui,vi, and the procedure returns a curve Nt associated to 
this vector field. 

> r2devlim3:=proc(M,p,q) 

> local prod, res , resl ,res2 , resul ,resu2 ,N,MMM: 

> MMM:=(r,s)->(op(M(r,s))) : 

> prod:=unapply([(eetal@(A(af@(ppsi3(3MMM@eeta3(3g@B(3ppsil))) (r,s)] ,t) : 

> res : =simplif y (subs (t=0,dif f (prod(t) , t) ) ) : 



> resl : =unapply (simplify (op (1 ,res) ) , (r , s) ) : 

> res2:=unapply(simplify(op(2,res)) , (r,s)) : 

> resul : =simplif y (mtaylor (resl (r , s) , [r=p, s=q] , 3) ) : 

> resu2 : =simplify (mtaylor (res2(r , s) , [r=p,s=q] ,3)) : 

> N:=unapply( [r+t*resul , s+t*resu2] , (r,s)) : 

> RETURN (N): 

> end: 

The procedure r2extract associates to a vector field given by a curve Nt{r2, S2) the components in dr^, ^2 dr2, 
S2dr2, r2S2dr2, ''I 1 '9r2 1 f^s2 1 '^2 c?s2 ; ^2 c?s2 ; ^252 5^2, ''i f^s2 ^ud 82832 in the Taylor expansion of the 
associated vector field at (0, 0). 

> r2extract : =proc (M) 

> local prod, res, resl, res2,w: 

> prod:=unapply(M(r,s) ,t) : 

> res : =simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply (op(l ,res) , (r , s) ) : 

> res2:=unapply(op(2,res) , (r,s)) : 

> w[l] :=simplify(resl(0,0)) : 

> w[2] :=simplify(subs({r=0,s=0},diff (resl(r,s) ,r))) : 

> w[3] :=simplify(subs({r=0,s=0},diff (resl(r,s) ,s))) : 

> w[4] :=simplify(subs({r=0,s=0},diff (resl(r,s) ,r,s))) : 

> w[5] :=simplify(subs({r=0,s=0},diff (resl(r,s) ,r$2))/2) : 

> w[6] :=simplify(subs({r=0,s=0},diff (resl(r,s) ,s$2))/2) : 

> w[7] :=simplify(res2(0,0)) : 

> w[8] :=simplify(subs({r=0,s=0},diff (res2(r,s) ,r))) : 

> w[9] :=simplify(subs({r=0,s=0},diff (res2(r,s) ,s))) : 

> w[10] :=simplify(subs({r=0,s=0},diff (res2(r,s) ,r,s))) : 

> w[ll] :=simplify(subs({r=0,s=0},diff (res2(r,s) ,r$2))/2) : 

> w[12] :=simplify(subs({r=0,s=0},diff (res2(r,s) ,s$2))/2) : 

> w: = [w[l] ,w[2] ,w[3] ,w[4] ,w[5] ,w[6] ,w[7] ,w[8] ,w[9] ,w[10] ,w[ll] ,w[12]] : 

> RETURN(convert(w,Vector)) : 

> end: 

The function pppsi-^ is the projection of the blow up of (0, 0)j^,z, (0, O)^^ and (0,0)^2,32; ^-^^d eeetai is its 
inverse. 

> pppsil : = (r , s)-> (r*s*s ,r*s*s"2) : 

> eeetal:=(y,z)->((y"2/z)/(z/y) , (z/y)) : 

The function pppsi2 is the projection of the blow up of (— ^, ^) , (0, j^^y ) and (0,0)^2,82; f-^id eeeta2 
is its inverse. 

> pppsi2:=(r,s)->(r*s*s-l/alpha, r*s*s*(s+l/(alpha-l))+l/alpha) : 

> eeeta2:=(y,z)->(((y+l/alpha)/((z-l/alpha)/(y+l/alpha)-l/(alpha-l)))/((z-l/alpha)/(y+l/alpha) 
-l/(alpha-l)) , ((z-l/alpha)/(y+l/alpha)-l/(alpha-l))) : 

The function pppsig is the projection of the blow up of ( — , — ) , (0, ) and (0,0)^2.^2' ^'^^ eeetas 
its inverse. 

> pppsi3:=(r,s)->(r*s*s+l/alpha, r*s*s*(s+l/(alpha+l))+l/alpha) : 

> eeeta3:=(y,z)->(((y-l/alpha)/((z-l/alpha)/(y-l/alpha)-l/(alpha+l)))/((z-l/alpha)/(y-l/alpha) 
-l/(alpha+l)) , ((z-l/alpha)/(y-l/alpha)-l/(alpha+l))) : 

We compute the points where we take the Taylor expansions. 

> limit (op (1 , [(eeeta2@(A@pppsil)) (r,s)] ) ,{r=l,s=0}) , limit (op (2, [(eeeta2@(A@pppsil)) (r ,s) 
]),{r=l,s=0}); 



IS 
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> limit(op(l , [(eeeta3@(A@f (3A)(3pppsil) (r ,s)] ) ,{r=l ,s=0}) , limit(op(2, [(eeeta3(3(A@f @A)(3pppsil) 
(r,s)]) ,{r=l,s=0}) ; 

(l + g)^ „ 

> limit (op (1 , [(eeetal@(A@f QAQf@A)(apppsil) (r,s)] ) ,{r=l,s=0}) , limit (op (2, [(eeetalQ(A(af OAQf OA) 
Opppsil) (r , s) ] ) , {r=l , s=0}) ; 

-1, 

The procedure r3devliml associates to a vector field Z on the blow up of at (0, 0)j^,z , {0,0)ui,vi and 
(0, 0)^2, S2 expressed in coordinates (ra, S3) the Taylor expansion of A^,Z of order 1 at a point (p, q). The result 



is viewed as a polynomial vector field in a neighborhood of {p, q) in the coordinates (ra, S3) which are this time 
the coordinates on blown up at (— -) , (0, and (0,0),.2^S2, and the procedure gives a curve 

Nt associated to this vector field. 

> r3devliml:=proc(M,p,q) 

> local prod, res ,resl ,res2 ,resul ,resu2 ,N,MMM: 

> MMM:=(r,s)->(op(M(r,s))) : 

> prod:=unapply([(eeeta2@(A@(pppsil(aMMM(aeeetal(3(B(3pppsi2)))) (r,s)] ,t) : 

> res :=simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply (simplify (op (1 ,res) ) , (r , s) ) : 

> res2:=unapply(simplify(op(2,res)) , (r,s)) : 

> resul : =simplif y (mtaylor (resKr , s) , [r=p, s=q] , 2) ) : 

> resu2 : =simplify (mtaylor (res2(r , s) , [r=p,s=q] ,2)) : 

> N: =unapply ( [r+t*resul , s+t*resu2] , (r , s) ) : 

> RETURN (N): 

> end: 

The procedure r3devlim2 associates to a vector field Z on the blow up of F'^ at (—-, -) , ( 0, -^-^ ) a-^d 

(0, 0)^2, S2 expressed in the coordinates (ra, S3) the Taylor expansion of {A(f))^, Z of order 1 at a point (p, q). The 
result is viewed as a polynomial vector field in a neighborhood of {p,q) in the coordinates (r3,S3) which are 

this time the coordinates of blown up at (i, , (O, ] and (0, 0)^2,821 and the procedure returns 

a curve Nt associated to this vector field. 

> r3devlim2:=proc(M,p,q) 

> local prod, res , resl ,res2 , resul ,resu2 ,N,MMM: 

> MMM: = (r,s)->(op(M(r,s))) : 

> prod:=unapply([(eeeta3@(A(Sf(a(pppsi2@MMM@eeeta2@g@B@pppsi3))) (r,s)] ,t) : 

> res : =simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply (simplify (op (1 ,res) ) , (r , s) ) : 

> res2:=unapply(simplify(op(2,res)) , (r,s)) : 

> resul : =simplify (mtaylor (resl (r , s) , [r=p, s=q] , 2) ) : 

> resu2 : =simplify (mtaylor (res2(r , s) , [r=p,s=q] ,2)) : 

> N : =unapply ( [r+t*resul , s+t*resu2] , (r , s) ) : 

> RETURN (N): 

> end: 

The procedure r3devlim3 associates to a vector field Z on the blow up of P^ at (-, -) , (0, ) and 

(0,0)^2,82 expressed in coordinates (tztS-^ the Taylor expansion of (A0)^ Z of order 1 at a point {p,q). The 
result is viewed as a polynomial vector field in a neighborhood of (p, q) in the coordinates {rs, S3) which are this 
time the coordinates on P^ blown up at (0, 0)y,z, (0, 0)ui,vi and (0, 0)^2,321 and the procedure returns a curve Nt 
associated to this vector field. 

> r3devlim3:=proc(M,p,q) 

> local prod, res, resl, res2, resul, resu2,N,MMM: 

> MMM:=(r,s)->(op(M(r,s))) : 

> prod:=unapply([(eeetal@(A(af@(pppsi3(3MMM@eeeta3(3g@B@pppsil))) (r,s)] ,t) : 

> res : =simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply (simplify (op (1 ,res) ) , (r , s) ) : 

> res2:=unapply(simplify(op(2,res)) , (r,s)) : 

> resul : =simplify (mtaylor (resl (r , s) , [r=p, s=q] , 2) ) : 

> resu2 : =simplify (mtaylor (res2(r , s) , [r=p,s=q] ,2)) : 

> N : =unapply ( [r+t*resul , s+t*resu2] , (r , s) ) : 

> RETURN (N): 

> end: 

The procedure r3extract associates to a vector field given by a curve Nt^r^^ss) the components of 9^3, 539^3 
and in the Taylor expansion of the associated vector field at (0,0). 

> r3extract : =proc (M) 

> local prod, res, resl, res2,w: 

> prod:=unapply(M(r,s) ,t) : 

> res : =simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply (op(l ,res) , (r , s) ) : 

> res2:=unapply(op(2,res) , (r,s)) : 

> w[l] :=simplify(resl(-l,0)) : 

> w[2] : =simplif y(subs({r=-l , s=0} ,dif f (resl (r , s) ,s))) : 

> w[3] :=simplify(res2(-l,0)) : 

> w: = [w[l] ,w[2] ,w[3]] : 

> RETURN(convert(w,Vector)) : 

> end: 



The function ppppsij^ is the projection of the blow up of (0,0)^^2, {0,0)ui,vi (OiO)r2,s2 s-^d (1,0)^3,83, and 
eeeetai is its inverse. 

> ppppsil:=(r,s)->((r*s+l)*s*s, (r*s+l) *s*s~2) : 

> eeeetal:=(y,z)->(((y-2/z)/(z/y)-l)/(z/y) , (z/y)) : 

The function ppppsi2 is the projection of the blow up at (—i, '(Oj^jTri) , (0, 0)r2,s2 a-^d ( — ^"^^i-* ,0) , 
and eeeeta2 is its inverse. 

> ppppsi2 : =(r , s)-> ( (r*s- (1/2) * (-1+alpha) "5/alpha"4) *s*s-l/alpha , (r*s- ( 1/2) * (-1+alpha) '5/ 
alpha'4) *s*s* (s+1/ (alpha- 1) ) +l/alpha) : 

> eeeeta2:=(y,z)->((((y+l/alpha)/((z-l/alpha)/(y+l/alpha)-l/(alpha-l)))/((z-l/alpha)/(y+l/ 
alpha)-l/(alpha-l) )+ (1/2)* (-1+alpha) -5/alpha~4)/( (z-l/alpha)/(y+l/alpha)-l/(alpha-l)) , ((z-1/ 
alpha)/(y+l/alpha)-l/(alpha-l))) : 

The function ppppsig is the projection of the blow up at ^) , (O, j^nrf) , (0:0)r2,s2 and ( ^^^^¥-,0), 
and eeeeta3 is its inverse. 

> ppppsi3:=(r,s)->((r*s+ (1/2)* (1+alpha) ~5/alpha"4)*s*s+l/alpha, (r*s+( 1/2) *( 1+alpha) "5/alpha~4) 
*s*s*(s+l/(alpha+l))+l/alpha) : 

> eeeeta3:=(y,z)->((((y-l/alpha)/((z-l/alpha)/(y-l/alpha)-l/(alpha+l)))/((z-l/alpha)/(y-l/ 
alpha)-l/(alpha+l))-(l/2)*(l+alpha)"5/alpha'4)/((z-l/alpha)/(y-l/alpha)-l/(alpha+l)) , ((z-1/ 
alpha) / (y- 1/alpha) - 1 / (alpha+1 ) ) ) : 

The function ppppsi^ is the projection of the blow up at (0,0)y.2, (0,0)ui,vi (0,0)^2.32 and (—1,0)^3,33, and 
eeeetai is its inverse. 

> ppppsi4:=(r,s)->((r*s-l)*s*s, (r*s-l) *s*s"2) : 

> eeeeta4:=(y,z)->(((y"2/z)/(z/y)+l)/(z/y) , (z/y)) : 

We compute the points where we take the Taylor expansions. 

> limit(op(l, [(eeeeta2@(A@ppppsil)) (r,s)] ) ,{r=0,s=0}) , limit(op(2, [(eeeeta2(S(A(Sppppsil) ) ( 
r,s)]),{r=0,s=0}); 

3{a ~ ~ + - a) 

> limit(op(l, [(eeeeta3@(A(Sf (SA)@ppppsil) (r,s)]) ,{r=0,s=0}) , limit(op(2, [(eeeeta3@(A@f @A)(Sppppsil) 
(r,s)]) ,{r=0,s=0}) ; 

-3(1 + a)^(~a - l + g^ + a^) 

1 ^ 



4a5 

> limit (op (1 , [(eeeeta4@(A@f OAOf @A)@ppppsil) (r, s)] ) , {r=0, s=0}) , limit (op (2, [(eeeeta4(a(A(af (SAOf OA) 
Oppppsil) (r,s)]) ,{r=0,s=0}) ; 

0, 

The procedure r4devliml associates to a vector field Z on the blow up of at {0,0)y^z , (0,0)„j_„j, (0,0)^2,82 
and (1, 0)^3, S3 expressed in the coordinates (r4, 34) the Taylor expansion of A^^Z of order at a point {p, q). The 
result is viewed as a polynomial vector field in a neighborhood of {p,q) in the coordinates (r^jSi) which are 

this time the coordinates on blown up at (— ^, ^) , (O, ;^^) , (0,0)^2,^2 and ( — ^°'2^J'J , ) , and 

the procedure returns a curve Nt associated to this vector field. 

> r4devliml:=proc(M,p,q) 

> local prod, res, resl,res2,resul,resu2,N,MMM: 

> MMM:=(r,s)->(op(M(r,s))) : 

> prod: =unapply( [(eeeeta2@(A(S(ppppsil(aMMM@eeeetal@(B@ppppsi2) )))(r,s)],t): 

> res : =simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) 

> resl : =unapply (simplify (op (1 ,res) ) , (r , s) ) 

> res2:=unapply(simplify(op(2,res)) , (r,s)) 

> resul : =simplif y (mtaylor (resl (r , s) , [r=p, s=q] , 1) ) : 

> resu2 : =simplify (mtaylor (res2(r , s) , [r=p, s=q] , 1) ) : 

> N : =unapply ( [r+t*resul , s+t*resu2] , (r , s) ) : 

> RETURN (N): 

> end: 

The procedure r4devlim2 associates to a vector field Z on the blow up of P^ at (— ^, ^) , (0, j^z^f ) i 
(0,0)r2,s2 and ( — ^"^'^^ ,0) expressed in the coordinates (r^^Si) the Taylor expansion of (Acj))^ Z of order 



at a point {p,q). The result is viewed as a polynomial vector field in a neighborhood of {p,q) in the coor- 



> r4devlim2:=proc(M,p,q) 

> local prod, res ,resl ,res2 ,resul ,resu2 ,N,MMM: 

> MMM:=(r,s)->(op(M(r,s))) : 

> prod:=unapply([(eeeeta3(3(A(af@(ppppsi2(aMMM@eeeeta2@g@B@ppppsi3))) (r,s)] ,t) : 

> res : =simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply (simplify (op (1 ,res) ) , (r , s) ) : 

> res2:=unapply(simplify(op(2,res)) , (r,s)) : 

> resul : =simplif y (mtaylor (resKr , s) , [r=p, s=q] , 1) ) : 

> resu2 : =simplify (mtaylor (res2(r , s) , [r=p, s=q] , 1) ) : 

> N:=unapply( [r+t*resul , s+t*resu2] , (r,s)) : 

> RETURN (N): 

> end: 

The procedure r4devlim3 associates to a vector field Z on the blow up of at (-, -) , (0, ^rr ) , 

(0,0)r2,s2 Slid ^ ^"^Ji , , and expressed in coordinates (r4,S4) the Taylor expansion of (Ac/))^ Z of order at 
a point {p, q) . The result is viewed as a polynomial vector field in a neighborhood of (p, q) in the coordinates 
(7-4,54) which are this time the coordinates on P^ blown up at (0, 0)y,z, (0,0)„j^^^, (0,0)^2,32 a-^d (—1,0)^3,33, 
and the procedure returns a curve Nt associated to this vector field. 

> r4devlim3:=proc(M,p,q) 

> local prod, res , resl ,res2 , resul ,resu2 ,N,MMM: 

> MMM:=(r,s)->(op(M(r,s))) : 

> prod:=unapply([(eeeeta4@(A@f<a(ppppsi3(3MMMQeeeeta3@g@B@ppppsi4))) (r,s)] ,t) : 

> res : =simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply (simplify (op (1 ,res) ) , (r , s) ) : 

> res2:=unapply(simplify(op(2,res)) , (r,s)) : 

> resul : =simplify (mtaylor (resl (r , s) , [r=p, s=q] , 1) ) : 

> resu2 : =simplify (mtaylor (res2(r , s) , [r=p, s=q] , 1) ) : 

> N : =unapply ( [r+t*resul , s+t*resu2] , (r , s) ) : 

> RETURN (N): 

> end: 

The procedure r4extract associates to a vector field given by a curve Nt{r4, S4) the components of dr^ and dg^ 

in the Taylor expansion of the vector field associated at (0, 0). 

> r4extract : =proc (M) 

> local prod, res, resl, res2,w: 

> prod:=unapply(M(r,s) ,t) : 

> res :=simplif y (subs (t=0 ,dif f (prod(t) , t) ) ) : 

> resl : =unapply (op(l ,res) , (r , s) ) : 

> res2:=unapply(op(2,res) , (r,s)) : 

> w[l] :=simplify(resl(0,0)) : 

> w[2] :=simplify(res2(0,0)) : 

> w: = [w[l] ,w[2]] : 

> RETURN(convert(w,Vector)) : 



6.3 Definition of the residues matrices 



Ri, denoted by in the article, is the (10 x 7)-matrix of dy, y dy, z dy, dz,y dz, zdz, dz in the basis (Ii)i<i<io- 

> R_l :=Transpose (Matrix ( [[-lambda [1] ,-1 ,2, 0,0,0, -lambda [4] "2 , -2*lambda [4] ,0,0] , [0,0,0,0,0, 
0,-3,0,0,0] , [0,0,0,0,0,0,0,0,-3,0] , [1 ,0,0, 0,0,0, lambda [4] "3,3*lambda[4] "2,2*lambda[5] "2,4* 
lambda [5]] , [0,0, -lambda [2] , -1 , -2*lambda [3] ,-2,-3*lambda[4] , -3 , -4*lambda [5] ,-4] , [0,0,0,0,0, 
0,2,0,0,0] , [0,0,0,0,0,0,0,0,-2,0]])); 




> end: 



R.l := 
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R2, denoted by ^ in the article, is the (10x7)-matrixof ydy, zdy, dz,ydz, zdz, dz in the basis (mi)i<i<io- 

> R_2 :=Transpose (Matrix ( [[-mu[l] ,-l,2,0,0,0,-mu[4] "2,-2*mu[4] ,0,0], [0,0, 0,0, 0,0, 3, 0,0,0], 
[0,0,0,0,0,0,0,0,3,0], [1 ,0,0, 0,0,0, -mu [4] ~3,-3*mu[4] ~2,2*mu[5] ~2,4*mu[5]] , [0,0,-mu[2] ,-1,-2* 
mu [3] , -2 , -3*mu [4] , -3 , -4*mu [5] , -4] , [0 , , , , , , -2 , , , 0] , [0 , , , , , , , , -2 , 0] ] ) ) ; 
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RRi is the (10 x 10)-matrix that gives the residues of d^^, Wii9„j, Vidm, Vidm, v^du^, dy-^, Uidy-^, 

vi dy^, vl dy^ in the basis ((i)i<i<io. 

> RR_l:=Matrix( [[0,0,0,0,0,0,0,0,0,0] ,[0,0,0,0,0,0,0,0,0,0], [1,0, 0,0, 0,0, -lambda [2] ,0,0,0] , 
[0,0,0,0,0,0,-1,0,0,0] , [0,0,l,0,0,0,-2*lambda[3] , 0, , 0] , [0 , 0, , ,0 , ,-2 , ,0 , 0] , [0 , -1 , , , 1 , 
0,-3*lambda[4] ,0,2,0] , [0,0,0,0,0,0,-3,0,0,0] , [0,0,0,-1 ,0, 1 ,-4*lambda[5] ,-2,0,2] , [0,0,0,0,0, 
0,-4,0,0,0]]) ; 































































1 

















-A2 





























-1 

















1 











-2 A3 





























-2 














-1 








1 





-3A4 





2 























-3 




















-1 





1 


-4 As 


-2 





2 




















-4 












RRRi is the (10 x 12)-niatrix that gives the residues of 9^2 1 ''2 dr2 1 ^2 1 ^2^2 9,-2 , f2^r2 ; ^2 9r2 ; , r2 , S2 ds,_ 
r2S2ds2, rjds2, ^2812 in the basis ([j)i<i<io- 

> RRR_l:=Matrix([[0,0,0,0,0,0,0,0,0,0,0,0] ,[0,0, 0,0, 0,0, 0,0, 0,0, 0,0], [0,0, 0,0, 0,0, 0,0, 0,0, 
0,0], [0,0, 0,0, 0,0, 0,0, 0,0, 0,0], [1 ,0,0, 0,0,0, -lambda [3] ,0,0, 0,0,0], [0,0, 0,0, 0,0, -1,0, 0,0,0, 



0] , [0,-1 , 1 ,0,0,0,-2*lambda[4] ,-1 , 1 ,0,0,0] , [0,0,0,0,0,0,-2,0,0,0,0,0] , [0,0,0,-1 , 1 , 1 , -3*lambda [5] , 
0,0,-1,1,1] , [0,0,0,0,0,0,-3,0,0,0,0,0]]); 



RRRJ 



' Order : 
> evalm(RRR_l) ; 



' 10 X 12 
'Data Type: 



(Matrix) 
' anything 



' Storage : 



' rectangular 
Fortran ^order 
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RRRRi is the (10 x 3)-niatrix that gives the residues of 9^3, ssdr^, ds^ in the basis ([i)i<i<io- 

> RRRR_l:=Matrix([ [0,0,0] , [0,0,0] , [0,0,0] , [0,0,0] , [0,0,0] , [0,0,0] , [1 ,0 , -lambda[4] ] , [0,0,-1] , 
[0,1, -2*lambda [5] ] , [0 , , -2] ] ) ; 



RRRR-l : — 
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RRRRRi is the (10 x 2)-matrix that gives the residues of dr^, dg^ in the basis ([i)i<i<io- 

> RRRRR_l:=Matrix([[0,0] , [0,0] , [0,0] , [0,0] , [0,0] , [0,0] , [0,0] , [0,0] , [1 , -lambda [5] ] , [0,-1]]) 



RRRRR^l '. — 
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-1 



P is the transition matrix from the algebraic basis (m.i)i<io to the geometric basis for P2. 



> P:=Matrix( [[-mu[l] ,1,0, 0,0, 0,0, 0,0,0], [-1,0, 0,0, 0,0, 0,0, 0,0], [2 , , 1 , -mu [2] ,0,0,0,0,0,0], 
[0,0,0,-1 ,0,0,0,0,0,0] , [0,0,0,-2*mu[3] ,l,-mu[3] ,0,0,0,0] , [0,0,0,-2,0,-1 ,0,0,0,0] , [-mu[4] "2, 
-mu[4] '3,0,-3*mu[4] ,0,-2*mu[4] ,1 ,-mu[4] ,0,0] , [-2*mu[4] ,-3*mu[4] "2,0,-3,0,-2,0,-1 ,0,0] , [0,2* 
mu[5] '■2,0,-4*mu[5] ,0,-3*mu[5] ,0,-2*mu[5] , 1 ,-mu[5]] , [0,4*mu[5] ,0,-4,0,-3,0,-2,0,-1]] ) ; 



P := 
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PP, denoted by in the article, is the transition matrix from the algebraic basis (li)i<i<io to the geometric 
basis (ei)i<i<io for Pi. 

> PP:=Matrix( [[-lambda [1] ,1,0,0,0,0,0,0,0,0] , [-1,0,0,0,0,0,0,0,0,0] , [2 , , 1 , -lambda [2] ,0,0, 
0,0,0,0] , [0,0,0,-1,0,0,0,0,0,0] , [0,0,0,-2*lambda[3] , 1 , -lambda [3] ,0,0,0,0] , [0,0,0,-2,0,-1 ,0, 

0. 0.01 , [-lambda [4] "2, lambda [4] "3,0,-3*lambda[4] ,0,-2*lambda[4] , 1, -lambda [4] ,0,0] , [-2*lambda[4] , 
3*lambda[4] "2,0,-3,0,-2,0,-1,0,0], [0,2*lambda[5] "2,0,-4*lambda[5] ,0,-3*lambda[5] ,0,-2*lambda[5] , 

1 , -lambda [5] ] , [0,4*lambda[5] , ,-4 , 0, -3 , 0, -2 , 0, -1] ] ) ; 



PP := 
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PPP is the transition matrix from the algebraic basis ((i)i<i<io to the geometric basis (ei)i<i<io for Pi. 
> PPP : =MatrixInverse (PP) : 



6.4 Construction of the matrix ^ 
6.4.1 Construction of a block matrix 

We input the coefficients of the matrix 

> col[l] :=[-mu[l] , -1, 2, 0, 0, 0, 2*mu[4]"2, 4*mu[4] , 0, 0, 0, 0, 0, 0, 0, 0, -4, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0]; 

coll := [-Ml, -1,2, 0,0, 0,2^1, 4^4, 0,0, 0,0, 0,0, 0,0, -4,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,-3, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,3,0] 

> col [2]: = 

> [1, 0, 0, 0, 0, 0, -mu[4]"3, -3*mu[4]"2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, -2]; 



coh := [1,0,0,0,0,0,-Ai!^, -3^^4,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,1,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-2] 

> col [3]: = 

> [0, 0, 1, 0, 0, 0, mu [4] "2, 2*mu[4], 0, 0, 0, 0, 0, 0, 0, 0, -2, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 1, 0]; 

cok [0, 0,1,0, 0, 0, nl, 2 /i4, 0, 0, 0, 0, 0, 0, 0, 0, -2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -2, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0] 

> col [4]: = 

> [0, 0, -mu[2], -1, 2*mu[3] "2-2*mu[3] , 4*mu[3]-2, mu[4], 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0]; 

cok := [0, 0, -/i2, -1, 2 - 2 /^3, 4/i3 - 2, ^4, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0,0] 

> col [5]: = 

> [0, 0, -2*mu[2], -2, (mu[3]-l)'2, 2*(mu[3]-l), 0, 0, 0, 0, 0, 0, -2, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0]; 

cok [0, 0, -2 Ai2, -2, (/ig - 1)^ 2 ^3 - 2, 0, 0, 0, 0, 0, 0, -2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0] 

> col [6]: = 

> [0, 0, 0, 0, mu[3] "2-iiiu[3] , 2*mu[3]-l, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0]; 

cok := [0,0,0,0,^^3 - M3,2/i3 - 1,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0] 

> col [7]: = 

> [0, 0, mu[2]'2, 2*mu[2], 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 1, 0, 0, 0] ; 

colj [0,0,/^^, 2 ^2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0] 

> col [8]: = 

> [0, 0, mu[2], 1, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0] ; 

cok ■■= [0,0,Ai2, 1,0, 0,0, 0,0, 0,0, 0,2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0 
,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

> col[9]:=[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 1, 0, 0, 0, 0] ; 

cok := [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,1,0, 0, 0, 0] 

> col [10]: = 

> [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0] ; 

coko [0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

The matrix PART is by definition if^-^ 

> PART: =MatrixMatrixMultiply (Transpose (MatrixC [col [1] ,col[2] ,col[3] ,col[4] ,col[5] ,col[6] , 
col [7] , col [8] , col [9] , col [10] ] ) ) , PPP) ; 

' 65 X 10 ' (Matrix) 

' Data Type : ' anything 

' Storage : ' rectangular 



PART 



' Order : ' Fortran _order 

> for k from 1 to 85 do 

> for j from 1 to 30 do 

> L[k][j]:=0 



> od 

> od: 



We construct the matrix ^ , whose upper right block =S is indeterminate. 

> for k from 1 to 10 do 

> for 1 from 1 to 10 do 

> L[k] [1+20] :=q[k,l] 

> od 

> od: 

> for k from 76 to 85 do 

> L[k] [k-65] :=1 

> od: 

> for k from 11 to 75 do 

> for j from 1 to 10 do 

> L[k] [j] :=PART[k-10, j] 

> od 

> od: 

> for k from 1 to 85 do 

> L[k] :=convert(L[k] ,list) 

> od: 

> MM:=Matrix( [seq(L[k] ,k=l. .85)]) : 
6.4.2 Computation of =S 

To compute the coefficients of =S, we consider the matrix R obtained by multiplying £2 at right by the transition 
matrix from the algebraic basis to the geometric basis of P2 • The matrix R is then the matrix 

:M^(P2,%eom) M^(A,^alg)- 

Computation of dy in the basis ((i)i<i<io- 

> y:='y': z:='z': 

> M: = (y,z)->[y+t,z] : 

> U:=MatrixMatrixMultiply(R_l,extract(devlim2(devlim2((devliml(M,-l/alpha,l/alpha)) , 1/alph 
1/alpha) ,0,0))) : 

> for k from 1 to 10 do 

> R[k] [1] :=U[k] ; 

> od: 

Computation of (^A(^A(^A)^ in the basis ([i)i<i<io- 

> y:='y': z:='z': 

> M: = (y,z)->[y,z+t] : 

> U:=MatrixMatrixMultiply(R_l,extract(devlim2(devlim2((devliml(M, -1/alpha, 1/alpha)) , 1/alph 
1/alpha) ,0,0))) : 

> for k from 1 to 10 do 

> R[k] [2] :=U[k] ; 

> od: 

Computation of {A4)A4)A)^, d^^ in the basis ([i)i<i<io. 

> u:='u': v:='v': 

> M:=(u,v)-> [u+t,v] : 

> U:=MatrixMatrixMultiply(RR_l,ulextract(uldevlim3(uldevlim2(uldevliml(M,0,l/(alpha-l)) , 0, 
l/(l+alpha)) ,0,0))) : 

> for k from 1 to 10 do 

> R[k] [3] :=U[k] ; 

> od: 

Computation of {A<pA<pA)^ d^i in the basis ([i)i<i<io- 

> u:='u': v:='v': 

> M:=(u,v)->[u,v+t] : 

> U: =MatrixMatrixMultiply(RR_l ,ulextract(uldevlim3(uldevlim2(uldevliml (M,0, l/(alpha-l) ) , 0, 
l/(l+alpha)) ,0,0))) : 

> for k from 1 to 10 do 

> R[k] [4] :=U[k] ; 

> od: 

Computation of 9^2 ™ the basis ([i)i<i<io. 

> r:='r': s:='s': 

> M:=(r,s)->[r+t,s] : 

> U : =MatrixMatrixMultiply (RRR_1 , r2extract (r2devlim3 (r2devlim2 (r2devliml (M, 0,0), 0,0), 0,0)) 
): 



> for k from 1 to 10 do 

> R[k] [5] :=U[k] ; 

> od: 

Computation of {A4>A(pA)^, ds^ in the basis ([i)i<i<io. 

> r:='r': s:='s': 

> M:=(r,s)->[r,s+t] : 

> U : =MatrixMatrixMultiply (RRR_1 , r2extract (r2devlim3 (r2devlim2 (r2devliml (M, 0,0), 0,0), 0,0)) 
): 

> for k from 1 to 10 do 

> R[k] [6] :=U[k] ; 

> od: 

Computation of {A(f)A(f)A)^ dr^ in the basis ([i)i<i<io. 

> r:='r': s:='s': 

> M: = (r,s)->[r+t,s] : 

> U : =MatrixMatrixMultiply (RRRR_ 1 , rSextract (rSdevlimS (r3devlim2 (rSdevliml (M , - ( 1/2) * (-1+alpha) 
~5/alpha"4,0) , (l/2)*(l+alpha) -5/alpha"4, 0) ,-1,0))) : 

> for k from 1 to 10 do 

> R[k] [7] :=U[k] ; 

> od: 

Computation of in the basis ((i)i<i<io- 

> r:='r': s:='s': 

> M:=(r,s)->[r,s+t] : 

> U:=MatrixMatrixMultiply(RRRR_l , rSextract (r3devlim3(r3devlim2 (rSdevliml (M,-( 1/2)* (-1+alpha) 
~5/alpha~4,0) , (l/2)*(l+alpha) "5/alpha"4, 0) ,-1,0))) : 

> for k from 1 to 10 do 

> R[k] [8] :=U[k] ; 

> od: 

Computation of {A(f)A(f)A)^, dc^ in the basis ([i)i<i<io- 

> r:='r': s:='s': 

> M:=(r,s)->[r+t,s] : 

> U : =MatrixMatrixMultiply (RRRRR_1 , r4extract (r4devlimS (r4devlim2 (r4devliml (M , (S/4) * (-1+alpha) 
"4* (alpha"S-alpha-alpha~2+l) /alpha's , 0) , - (S/4) * (1+alpha) "4* (-alpha-l+alpha"S+alpha~2) /alpha~5 , 
0) ,0,0))) : 

> for k from 1 to 10 do 

> R[k] [9] :=U[k] ; 

> od: 

Computation of {A(f)A4iA)^ dd^ in the basis ([i)i<i<io- 

> r: = 'r' : s: = 's' : 

> M:=(r,s)->[r,s+t] : 

> U : =MatrixMatrixMultiply (RRRRR_ 1 , r4extract (r4devlimS (r4devlim2 (r4devliml (M , (S/4) * (-1+alpha) 
"4* (alpha"S-alpha-alpha'2+l) /alpha'5 , 0) , - (S/4) * (1+alpha) "4* (-alpha-l+alpha"S+alpha'2) /alpha'5 , 
0),0,0))): 

> for k from 1 to 10 do 

> R[k] [10] :=U[k] ; 

> od: 

We construct the matrix R, then we compute =S. 

> for k from 1 to 10 do 

> R [k] : =convert (R [k] , list) 

> od: 

> TEMP:=Matrix( [seq(R[k] ,k=l. .10)]): 

> Q:=MatrixMatrixMultiply(TEMP, MatrixInverse(P)) : 

We can display the coefficients of ^. 

> simplify(q[9,4]) ; 



4(23 + a-6A5)/X2 



6.5 The 8 vectors of V{Di) that form the columns of the matrix ^ 



6.5.1 Columns of (eight (10 x l)-vectors) 

ai 

> prol:=[-lambda[l] ,-1 ,2,0,0,0,-lambda[4] "2,-2*lambda[4] ,0,0] ; 

> a [1] : =convert (prol , Vector) : 

prol := [-Al,-l,2,0,0,0,-A4^-2A4,0,0] 

> pro2: = [0,0,0,0,0,0,-3,0,0,0] ; 

> a[2] :=convert(pro2, Vector) : 

pro2 ~ [0,0,0,0,0,0,-3,0,0,0] 
as 

> pro3: = [0,0,0,0,0,0,0,0,-3,0] ; 

> a [3] : =convert (pro3 , Vector) : 

proS := [0,0,0,0,0,0,0,0,-3,0] 

> pro4:=[l ,0,0, 0,0,0, lambda [4] ~3,3*lambda[4] "2 , 2*lambda[5] "2,4*lambda[5] ] ; 

> a[4] :=convert(pro4,Vector) : 

pro4 [1,0,0,0,0,0,A4',3A4',2A5',4A5] 
as 

> pro5 : = [0 , , -lambda [2] , -1 , -2*lambda [3] , -2 , -3*lambda [4] , -3 , -4*lambda [5] , -4] ; 

> a [5] : =convert (pro5 , Vector) : 

pro5 := [0,0,-A2,-l,-2A3,-2,-3A4,-3,-4A5,-4] 
ae 

> pro6: = [0,0,0,0,0,0,2,0,0,0] ; 

> a [6] : =convert (pro6 , Vector) : 

pro6 := [0,0,0,0,0,0,2,0,0,0] 
a? 

> pro7: = [0,0,0,0,0,0,0,0,0,0] ; 

> a [7] : =convert (pro7 , Vector) : 

pro? := [0,0,0,0,0,0,0,0,0,0] 
as 

> pro8: = [0,0,0,0,0,0,0,0,0,0] ; 

> a [8] : =convert (pro8 , Vector) : 

pro8 ~ [0,0,0,0,0,0,0,0,0,0] 

6.5.2 Columns of (eight (10 x l)-vectors) 

ai 

> M: = (y,z)->(t+y,z) : 

> aa[l] :=MatrixMatrixMultiply(R_2,champ4(M)) : 

a2 

> M:=(y,z)->((l+t)*y,z) : 

> aa[2] :=MatrixMatrixMultiply(R_2,champ4(M)) : 

as 

> M: = (y,z)->(y+t*z,z) : 

> aa[3] :=MatrixMatrixMultiply(R_2,champ4(M)) : 



a4 

> M: = (y,z)->(y,t+z) : 

> aa [4] : =Matr ixMatr ixMultiply (R_2 , champ4 (M) ) : 

> M:=(y,z)->(y,t*y+z) : 

> aa [5] : =MatrixMatr ixMultiply (R_2 , champ4 (M) ) : 

ae 

> M:=(y,z)->(y, (l+t)*z) : 

> aa [6] : =MatrixMatr ixMultiply (R_2 , champ4 (M) ) : 

ar 

> M:=(y,z)->(y/(l-t*y) ,z/(l-t*y)) : 

> aa [7] : =MatrixMatr ixMultiply (R_2 , champ4 (M) ) : 

OS 

> M:=(y,z)->(y/(l-t*z) ,z/(l-t*z)) : 

> aa [8] : =MatrixMatr ixMultiply (R_2 , champ4 (M) ) : 

6.5.3 Columns of (eight (10 x l)-vectors) 

ai 

> M: = (y,z)->(t+y,z) : 

> aaa[l] :=MatrixMatrixMultiply(R_2,chainp3(M)) : 

02 

> M:=(y,z)->((l+t)*y,z) : 

> aaa[2] :=MatrixMatrixMultiply(R_2,champ3(M)) : 

as 

> M: = (y,z)->(y+t*z,z) : 

> aaa[3] :=MatrixMatrixMultiply(R_2,champ3(M)) : 

> M: = (y,z)->(y,t+z) : 

> aaa[4] :=MatrixMatrixMultiply(R_2,champ3(M)) : 

05 

> M: = (y,z)->(y,t*y+z) : 

> aaa[5] :=MatrixMatrixMultiply(R_2,champ3(M)) : 

> M:=(y,z)->(y, (l+t)*z) : 

> aaa[6] :=MatrixMatrixMultiply(R_2,champ3(M)) : 

a? 

> M:=(y,z)->(y/(l-t*y) ,z/(l-t*y)) : 

> aaa[7] :=MatrixMatrixMultiply(R_2,chainp3(M)) : 

as 

> M:=(y,z)->(y/(l-t*z) ,z/(l-t*z)) : 

> aaa[8] :=MatrixMatrixMultiply(R_2,champ3(M)) : 



6.5.4 Construction of the matrix 



> mise:=proc(k: : integer) 

> local u,v,w,h: 

> u:=convert (Transpose (a [k] ) ,list) : 

> v:=convert (Transpose (aa[k] ) ,list) : 

> w : =convert (Transpose (aaa [k] ) , list) : 

> h:=[op(u) ,op(v) ,op(w)] : 

> RETURN ( convert (h, Vector)) 

> end: 

> for k from 1 to 8 do 

> VectDef [k] :=iiiise(k) 

> od: 

6.6 The 63 vectors of V{D2) that form the columns of the matrix ^ 

6.6.1 Vectors of type a (we compute the columns of the matrices ^/fl et JK^') 

Columns of ,yK^ {eight (10 x l)-vectors) 

Aai 

> M:=(y,z)->(t+y,z) : 

> aA[l] :=MatrixMatrixMultiply(R_l,champl(M)) : 

Aa2 

> M:=(y,z)->((l+t)*y,z) : 

> aA[2] :=MatrixMatrixMultiply(R_l,champl(M)) : 

> M: = (y,z)->(y+t*z,z) : 

> aA[3] :=MatrixMatrixMultiply(R_l,champl(M)) : 

Aa4 

> M: = (y,z)->(y,t+z) : 

> aA[4] :=MatrixMatrixMultiply(R_l,champl(M)) : 

Aa^ 

> M: = (y,z)->(y,t*y+z) : 

> aA[5] : =MatrixMatrixMultiply(R_l , champl(M) ) : 

Aae 

> M:=(y,z)->(y, (l+t)*z) : 

> aA[6] :=MatrixMatrixMultiply(R_l,champl(M)) : 

Aar 

> M:=(y,z)->(y/(l-t*y) ,z/(l-t*y)) : 

> aA[7] :=MatrixMatrixMultiply(R_l,champl(M)) : 

Aas 

> M:=(y,z)->(y/(l-t*z) ,z/(l-t*z)) : 

> aA[8] :=MatrixMatrixMultiply(R_l,champl(M)) : 

Columns of .yVg, to which we add 55 zeros 
Aai 

> pro 1 : = [-mu [1], -1,2, 0,0,0, -mu [4] " 2 , -2*mu [4], 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(prol) ; 

> aaA [1] : =convert (prol , Vector) : 

prol := [-/^i, -1,2, 0,0, 0,-^4, -2 ^4, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 

Aa2 



> pro2 := [0,0, 0,0, 0,0, 3, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro2) ; 

> aaA [2] : =convert (pro2 , Vector) : 

pro2 := [0,0,0,0,0,0,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 
,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 

> pro3:= [0,0, 0,0, 0,0, 0,0, 3, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro3) ; 

> aaA [3] : =convert (pro3 , Vector) : 

pro3 := [0,0,0,0,0,0,0,0,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0 
,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 

Aa4 

> pro4:=[l,0,0,0,0,0,-mu[4] "3,-3*mu[4] "2,2*mu[5] ~2,4*mu[5] ,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro4) ; 

> aaA [4] : =convert (pro4 , Vector) : 

pro4 [l,0,0,0,0,0,-/i|,-3Afi2/i^,4^5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 

Aa5 

> pro5 : = [0 , , -mu [2] , -1 , -2*mu [3] , -2 , -3*mu [4] , -3 , -4*mu [5], -4, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro5) ; 

> aaA [5] : =convert (pro5 , Vector) : 

pro5 [0,0,-Ai2, -1,-2/13,-2, -3/^4, -3, -4/^5, -4, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
6,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 

Aae 

> pro6 :=[0, 0,0, 0,0, 0,-2, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro6) ; 

> aaA [6] : =convert (pro6 , Vector) : 

pro6 := [0, 0, 0, 0, 0, 0, -2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
Aar 

> pro7:= [0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro7) ; 

> aaA [7] : =convert (pro7 , Vector) : 

pro7 := [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 

Aas 

> pro8 := [0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro8) ; 

> aaA [8] : =convert (pro7 , Vector) : 



pro8 ~ [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 



65 

Columns of ,yVJ' {eight (10 x l)-vectors) 
Aai 

> M: = (y,z)->(t+y,z) : 

> aaaA[l] :=simplify(MatrixMatrixMultiply(R_2,champ2(M))) : 

Aa2 

> M:=(y,z)->((l+t)*y,z) : 

> aaaA[2] :=simplify(MatrixMatrixMultiply(R_2,champ2(M))) : 

> M: = (y,z)->(y+t*z,z) : 

> aaaA[3] :=simplify(MatrixMatrixMultiply(R_2,champ2(M))) ; 

Aa4 

> M: = (y,z)->(y,t+z) : 

> aaaA[4] :=simplify(MatrixMatrixMultiply(R_2,champ2(M))) : 

Aa5 

> M:=(y,z)->(y,t*y+z) : 

> aaaA[5] :=simplify(MatrixMatrixMultiply(R_2,champ2(M))) : 

Aae 

> M:=(y,z)->(y, (l+t)*z) : 

> aaaA[6] :=simplify(MatrixMatrixMultiply(R_2,champ2(M))) : 

Aar 

> M:=(y,z)->(y/(l-t*y) ,z/(l-t*y)) : 

> aaaA[7] :=simplify(MatrixMatrixMultiply(R_2,champ2(M))) : 

Aas 

> M:=(y,z)->(y/(l-t*z) ,z/(l-t*z)) : 

> aaaA[8] :=simplify(MatrixMatrixMultiply(R_2,champ2(M))) : 

Glueing 

We form the first eight vectors of type a of y (2)2) ■ 

> for k from 1 to 8 do 

> for j from 1 to 10 do 

> VectDef_a[k] [j] :=aA[k] [j] 

> od 

> od: 

> for k from 1 to 8 do 

> for j from 11 to 75 do 

> VectDef_a[k] [j] :=aaA[k] [j-10] 

> od 

> od: 

> for k from 1 to 8 do 

> for j from 76 to 85 do 

> VectDef_a[k] [j] :=aaaA[k] [j-75] 

> od 

> od: 

> for k from 1 to 8 do 

> VectDef_a[k] :=convert(VectDef_a[k] ,list) 

> od: 

> for k from 1 to 8 do 

> VectDef_a[k] : =convert (VectDef _a [k] .Vector) 

> od: 



6.6.2 Vectors of type b (we compute the columns of the matrices yi^, ,Ah and .yV^') 

Columns of ^A^ {twenty Eve (10 x l)-vectors) 

> Mg:=(y,z,p,q)->(y+t*y"p/z"q,z) ; 

tyP 

Mg := {y,z,p,q) + —,z 

zi 

> M:=(y,z)->Mg(y,z,p,q) ; 

M := {y,z) i-^ Mg{y,z,p,q) 



> bA:=unapply(MatrixMatrixMultiply(R_l, champ Km)) , (p,q)) : 

Colums of o/fb to which we add a (25 x 25) -identity block then 30 zeros 
Abo,i 

> prol:=[0,0,0,0,0,0,2*mu[4] "5, 10*mu[4] "4,3*mu[5] ~3,9*mu[5] '2, 1,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(prol) ; 

> bbA[0, 1] :=convert(prol , Vector) : 

> bbA[0,l] [11] ; 

prol := [0, 0, 0, 0, 0, 0, 2 10 3 //f, 9 /x^, 1,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> pro2 : = [-1 , , , , , , mu [4] "3 , 3*mu [4] "2 , -2*mu [5] "2 , -4*mu [5], 0,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro2) ; 

> bbA[l , 1] :=convert(pro2, Vector) : 

> bbA[l,l] [12] ; 

pro2 [-1, 0,0,0,0,0,/ii 3 /Lt4, -2/^^,-4/^5, 0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

Ab2,l 

> pro3 : = [0 , , 2*mu [2] , 2 , 3*mu [3] , 3 , 4*mu [4] , 4 , 5*mu [5], 5, 0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro3) ; 

> bbA [2,1] : =convert (pro3 , Vector) : 

> bbA [2,1] [13] ; 

pro3 := [0, 0, 2 /i2, 2, 3 /i3, 3, 4 /i4, 4, 5 /is, 5, 0,0, 1,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
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Abo,2 

> pro4 := [0,0, 0,0, 0,0, -9*mu [4] "8 , -72*mu [4] "7, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro4) ; 

> bbA [0,2] :=convert(pro4, Vector) : 

> bbA [0,2] [14] ; 

pro4 [0,0, 0,0, 0,0, -9/4^,-72/*;^, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 



> pro5 : = [0,0, 0,0, 0,0, -3*mu [4] "6 , -18*mu [4] "5,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro5) ; 

> bbA [1,2] : =convert (pro5 , Vector) : 

> bbA [1,2] [15] ; 

pro5 [0,0,0,0,0,0,-3/1^,-18/^^0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
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> pro6:=[-l/mu[l] , l/mu[l] '2,0,0,0,0,-mu[4] "4,-4*mu[4] "3,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro6) ; 

> bbA [2, 2] : =convert (pro6 , Vector) : 

> bbA [2, 2] [16] ; 

pro6 := [-/t];\/t];^0,0,0,0,-/i|,-4/i|,0,0,0,0,0,0,0,l,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
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> pro7 := [0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro7) ; 

> bbA [3, 2] : =convert (pro7 , Vector) : 

> bbA [3, 2] [17] ; 

pro! := [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 

65 
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> pro8:=[0,0,0,0,0,0,52*mu[4] "11 , 572*mu [4] "10,-28*mu[5] "6 , -168*mu [5] "5,0,0,0,0,0,0,0,1,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0]; 

> nops(pro8) ; 

> bbA [0,3] : =convert (pro8 , Vector) : 

> bbA [0,3] [18] ; 

pro8 := [0, 0, 0, 0, 0, 0, 52 /tl\ 572 /^f , -28 /t^, -168 /i^, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> pro9:=[0,0,0,0,0,0, 15*mu[4] "9, 135*mu[4] "8, 12*mu[5] "5,60*mu[5] "4,0,0,0,0,0,0,0,0,1,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(prol) ; 

> bbA [1,3] : =convert (pro9 , Vector) : 

> bbA [1,3] [19] ; 

pro9 := [0, 0, 0, 0, 0, 0, 15 /i|, 135 fil, 12 /i^, 60 /i^ 0, 0, 0, 0, 0, 0, 0, 0,1,0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,6,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 



Ab2,3 

> prolO: = [0,0,0,0,0,0,4*mu[4] '■7,28*mu[4] ~6,-5*mu[5] "4,-20*mu[5] "3,0,0,0,0,0,0,0,0,0,1,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(prolO) ; 

> bbA[2,3] :=convert (pro 10, Vector) : 

> bbA[2,3] [20] ; 

prolO := [0, 0, 0, 0, 0, 0, 4/^!^, 28 -5 Afg, -20 /i^, 0,0, 0,0, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

Ab3,3 

> proll : = [-l/iiiu[l] '2,2/mu[l] "3,0,0,0,0,mu[4] "5,5*mu[4] "4,2*niu[5] "3,6*mu[5] "2,0,0,0,0,0,0, 
0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0] ; 

> nops(proll) ; 

> bbA[3,3] :=convert(proll, Vector) : 

> bbA[3,3] [21] ; 

proll [-^^^ 2 0,0, 0,0, /it, 5 Ait 2^5, 6^5, 0,0, 0,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 

0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> pro 12 := [0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(prol2) ; 

> bbA[4,3] :=convert(prol2, Vector) : 

> bbA[4,3] [22] ; 

prol2 := [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 


,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

Abo,i 

> prol3 := [0,0, 0,0, 0,0, -340*mu [4] " 14 , -4760*mu [4] "13,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(prol3) ; 

> bbA[0,4] :=convert(prol3, Vector) : 

> bbA[0,4] [23] ; 

prol3 := [0, 0, 0, 0, 0, 0, -340 nl'^, -4760 /if, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

AbiA 

> prol4 := [0,0, 0,0, 0,0, -91*mu [4] " 12 , -1092*mu [4] "11,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(prol4) ; 

> bbA[l,4] :=convert (pro 14, Vector) : 

> bbA[l,4] [24] ; 

prol4 [0,0,0,0,0,0,-91/if ,-1092 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 



Ab2,4 

> prol5:=[0,0,0,0,0,0,-22*mu[4] "10,-220*mu[4] "9,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(prol5) ; 

> bbA[2,4] :=convert(prol5, Vector) : 

> bbA[2,4] [25] ; 

prol5 [0, 0, 0, 0, 0, 0, -22 /x^", -220 nl, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> prol6 := [0,0, 0,0, 0,0, -5*mu [4] "8 , -40*mu [4] "7, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,1, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(prol6) ; 

> bbA[3,4] :=convert (pro 16, Vector) : 

> bbA[3,4] [26] ; 

prol6 := [0, 0, 0, 0, 0, 0, -5 nl, -40 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,1,0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> prol7:=[-l/mu[l] "3,3/mu[l] "4,0,0,0,0,-mu[4] "6,-6*mu[4] "5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(prol7) ; 

> bbA[4,4] :=convert (pro 17, Vector) : 

> bbA[4,4] [27] ; 

prol7 ~ [-^ir^, 3 A^r^, 0,0, 0,0,-^^^,-6^^4,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> pro 18 := [0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(prol8) ; 

> bbA[5,4] :=convert (pro 18, Vector) : 

> bbA[5,4] [28] ; 

prol8 ~ [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 

65 
1 

> prol9:=[0,0,0,0,0,0,2394*mu[4] "17,40698*mu[4] "16,429*mu[5] "9 , 3861*mu [5] "8,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0] ; 

> nops(prol9) ; 

> bbA[0,5] :=convert(prol9, Vector) : 

> bbA[0,5] [29] ; 

prol9 := [0, 0, 0, 0, 0, 0, 2394 /if , 40698 /xf , 429 /i|, 3861 /x^, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,0,1, 0, 0, 0, 0, 
0,0,0,0,0,0,0,6,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 



> pro20:=[0,0,0,0,0,0,612*mu[4] *15,9180*mu[4] "14,-165*mu[5] "8 , -1320*mu [5] "7,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0] ; 

> nops(pro20) ; 

> bbA [1,5]: =convert (pro20 , Vector) : 

> bbA [1,5] [30] ; 

pro20 := [0, 0, 0, 0, 0, 0, 612 /if , 9180 A^l'', -165 ^if, -1320 /i^, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

-462,5 

> pro21 := [0,0, 0,0, 0,0, 140*mu [4] " 13 , 1820*mu [4] ' 12 , 60*mu [5] "7 , 420*mu [5] "6, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0]; 

> nops(pro21) ; 

> bbA [2,5] : =convert (pro21 , Vector) : 

> bbA [2, 5] [31] ; 

pro21 [0, 0, 0, 0, 0, 0, 140 , 1820 /if , 60 /i^, 420 /if, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,1,0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

^63,5 

> pro22:=[0,0,0,0,0,0,30*mu[4] "ll,330*mu[4] "10,-21*mu[5] "6,-126*mu[5] "5,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0]; 

> nops(pro22) ; 

> bbA [3,5] : =convert (pro22 , Vector) : 

> bbA [3 , 5] [32] ; 

pro22 := [0, 0, 0, 0, 0, 0, 30 330 /Ltf , -21 /i|, -126 /i|, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

^64,5 

> pro23:=[0,0,0,0,0,0,6*mu[4] "9,54*mu[4] '8,7*mu[5] '5,35*mu[5] "4,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro23) ; 

> bbA [4,5] : =convert (pro23 , Vector) : 

> bbA [4, 5] [33] ; 

pro23 := [0, 0, 0, 0, 0, 0, 6 /i^^, 54 /f^, 7 /i^, 35 /i^, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

^65,5 

> pro24:=[-l/mu[l] "4,4/mu[l] "5,0,0,0,0,mu[4] '7,7*mu[4] "6,-2*mu[5] "4,-8*mu[5] "3,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0] ; 

> nops(pro24) ; 

> bbA [5,5] : =convert (pro24 , Vector) : 

> bbA [5 , 5] [34] ; 

pro24 [-^r^, 4/iJ;^ 0, 0, 0, 0, fil, 7 fil, -2 fi^ -8 /i^, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,1,0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 



65 



1 



> pro25:= [0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro25) ; 

> bbA [6,5] : =convert (pro25 , Vector) : 

> bbA [6, 5] [35] ; 

pro25 := [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 

65 
1 

Columns of JV^' {twenty Eve (10 x l)-vectors) 

> Mg:=(y,z,p,q)->(y+t*y'p/z"q,z) ; 

tyP 

Mg := {y,z,p,q) y-^ y + —,z 



> M:=(y,z)->Mg(y,z,p,q) ; 

M ~ {y,z) Mg{y,z,p,q) 

> bbbA:=unapply(MatrixMatrixMultiply(R_2,champ2(M)) , (p,q)) : 
Glueing 

We form the 25 vectors of type b of V{'D2)- 

> for q from 1 to 5 do 

> for p from to q+1 do 

> for j from 1 to 10 do 

> u[p,q] [j] :=bA(p,q) [j] 

> od 

> od 

> od: 

> for q from 1 to 5 do 

> for p from to q+1 do 

> for j from 11 to 75 do 

> u[p,q] [j] :=bbA[p,q] [j-10] 

> od 

> od 

> od: 

> for q from 1 to 5 do 

> for p from to q+1 do 

> for j from 76 to 85 do 

> u [p , q] [ j ] : =bbbA (p , q) [ j -75] 

> od 

> od 

> od: 

> for q from 1 to 5 do 

> for p from to q+1 do 

> u[p,q] :=convert(u[p,q] ,list) 

> od 

> od: 

> q:=l: 

> for p from to q+1 do 

> VectDef _b [p+1] : =u [p , q] 

> od: 

> q:=2: 

> for p from to q+1 do 

> VectDef _b [p+4] : =u [p , q] 

> od: 

> q:=3: 

> for p from to q+1 do 

> VectDef _b [p+8] : =u [p , q] 

> od: 

> q:=4: 

> for p from to q+1 do 

> VectDef _b [p+13] : =u [p , q] 



> od: 

> q:=5: 

> for p from to q+1 do 

> VectDef _b [p+19] : =u [p , q] 

> od: 

> p:='p': 

> q:='q': 



6.6.3 Vectors of type c (we compute the columns of the matrices J^J, ^ and ^") 

Columns of {twenty five (10 x l)-vectors) 

> Mg:=(y,z,p,q)->(y,z+t*y"p/z"q) ; 

Mg := {y,z,p,q) ^y,z+ — 

> M:=(y,z)->Mg(y,z,p,q) ; 

M := {y,z) i-^ Mg{y,z,p,q) 

> cA:=unapply(MatrixMatrixMultiply(R_l,champl(M)) , (p,q)) : 

Columns of ^/(^ to which we add 25 zeros, a (25 x 25)-identity block then 5 zeros 

> pro 1: = [0,0, 0,0, 0,0, 4*mu [4] " 6 , 24*mu [4] "5, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(prol) ; 

> ccA [0,1] : =convert (prol , Vector) : 

> ccA[0,l] [36] ; 

prol := [0, 0, 0,0,0,0,4/i^, 24 /x^, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,1,0,0,0,0,0,0,0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 

65 
1 

> pro2:=[l/mu[l] ,-l/mu[l] '2,0,0,0,0,mu[4] ~4,4*mu[4] "3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro2) ; 

> ccA [1,1]: =convert (pro2 , Vector) : 

> ccA[l,l] [37] ; 

pro2 := [//^^-^^^0,0,0,0,/xt4/i|, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 1,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 

1 

AC2,1 

> pro3:= [0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 1,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro3) ; 

> ccA [2,1] : =convert (pro3 , Vector) : 

> ccA[2,l] [38] ; 

proS := [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> pro4:=[0,0,0,0,0,0,-25*mu[4] '9,-225*mu[4] "8,-18*mu[5] "5,-90*mu[5] "4,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 



> nops(pro4) ; 

> ccA[0,2] :=convert(pro4, Vector) : 

> ccA[0,2] [39] ; 

pro4 := [0, 0, 0, 0, 0, 0, -25 -225 nl, -18 nl, -90 nt, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> pro5:=[0,0,0,0,0,0,-5*mu[4] "7,-35*mu[4] "6,6*mu[5] "4,24*mu[5] "3,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro5) ; 

> ccA [1,2] : =convert (pro5 , Vector) : 

> ccA[l,2] [40] ; 

pro5 := [0,0,0,0,0,0,-5^*1,-35^^,6/1^,24/^^,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> pro6:=[l/mu[l] "2,-2/mu[l] "3,0,0,0,0,-mu[4] "5,-5*mu[4] "4,-2*mu[5] "3,-6*mu[5] "2,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0] ; 

> nops(pro6) ; 

> ccA[2,2] : =convert (pro6 , Vector) : 

> ccA[2,2] [41] ; 

pro6 ~ [/<^^-2/i^^0,0,0,0,-//|,-5//t-2/ii,-6/i^,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,1,0,6,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> pro7: = [0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro7) ; 

> ccA[3,2] : =convert (pro7 , Vector) : 

> ccA[3,2] [42] ; 

pro? := [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 

65 
1 

> pro8:=[0,0,0,0,0,0, 182*mu[4] "12,2184*mu[4] "11,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro8) ; 

> ccA[0,3] : =convert (pro8 , Vector) : 

> ccA[0,3] [43] ; 

pro8 ~ [0,0,0,0,0,0,182/if ,2184 /i^\ 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> pro9 := [0,0, 0,0, 0,0, 33*mu [4] " 10 , 330*mu [4] "9, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 



> nops(prol) ; 

> ccA [1,3] : =convert (pro9 .Vector) : 

> ccA[l,3] [44] ; 

pro9 := [0, 0, 0, 0, 0, 0, 33 nl° , 330 fil, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> pro 10 : = [0 , , , , , , 6*mu [4] "8 , 48*mu [4] "7, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(prolO) ; 

> ccA[2,3] :=convert(prolO, Vector) : 

> ccA[2,3] [45] ; 

prolO := [0,0,0,0,0,0,6/1^,48^*1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 

65 
1 

> proll:=[l/mu[l] "3,-3/mu[l] "4,0,0,0,0,mu[4] "6,6*mu[4] "5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(proll) ; 

> ccA[3,3] :=convert(proll, Vector) : 

> ccA[3,3] [46] ; 

proll := [/lJ'^ -3 /i]"'',0,0,0,0,/z2, 6/^1,0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> pro 12:= [0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(prol2) ; 

> ccA[4,3] :=convert(prol2, Vector) : 

> ccA[4,3] [47] ; 

prol2 := [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> prol3:=[0,0,0,0,0,0,-1428*mu[4] "15,-21420*mu[4] "14,330*mu[5] '8,2640*mu[5] "7, 

> 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0]; 

> nops(prol3) ; 

> ccA[0,4] :=convert(prol3, Vector) : 

> ccA[0,4] [48] ; 

prolS [0, 0, 0, 0, 0, 0,-1428 fiY', -21420 330 /i^, 2640 14, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> prol4:=[0,0,0,0,0,0,-245*mu[4] "13,-3185*mu[4] "12,-96*mu[5] "7 , -672*mu [5] "6, 

> 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0]; 



> nops(prol4) ; 

> ccA[l ,4] :=convert (pro 14, Vector) : 

> ccA[l,4] [49] ; 

prol4 := [0, 0, 0, 0, 0, 0, -245 fif, -3185 , -96 nl, -672 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> pro 15 : = [0 , , , , , , -42*mu [4] " 1 1 , -462*mu [4] " 10 , 28*mu [5] "6 , 168*mu [5] "5, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0]; 

> nops(prol5) ; 

> ccA[2,4] :=convert(prol5, Vector) : 

> ccA[2,4] [50] ; 

prol5 := [0, 0, 0, 0, 0, 0, -42 -462 nl^, 28 168 ^^1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> prol6:=[0,0,0,0,0,0,-7*mu[4] ~9,-63*mu[4] ~8,-8*mu[5] ~5,-40*mu[5] "4,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(prol6) ; 

> ccA[3,4] :=convert(prol6, Vector) : 

> ccA[3,4] [51] ; 

prol6 := [0, 0, 0, 0, 0, 0, -7yLt^, -63^*4, -8 /Lt^, -40 Ai^, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> prol7:=[l/mu[l] "4,-4/mu[l] ~5,0,0,0,0,-mu[4] "7,-7*mu[4] "6,2*mu[5] "4,8*mu[5] "3,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0, 
0,0,0,0,0] ; 

> nops(prol7) ; 

> ccA[4,4] :=convert(prol7, Vector) : 

> ccA[4,4] [52] ; 

prol7 := [Ai^'',-4^J;^0,0,0,0,-/i|, -7/1^,2/^^, 8 ^^,0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

^C5,4 

> pro 18 := [0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(prol8) ; 

> ccA[5,4] :=convert(prol8, Vector) : 

> ccA[5,4] [53] ; 

prol8 := [0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 



> prol9 := [0,0, 0,0, 0,0, 1 1704*mu [4] ' 18 , 210672*mu [4] "17, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0]; 

> nops(prol9) ; 

> ccA[0,5] :=convert(prol9, Vector) : 

> ccA[0,5] [54] ; 

prol9 ~ [0, 0, 0, 0, 0, 0, 11704 /if , 210672 /if, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> pro20 := [0,0, 0,0, 0,0, 1938*mu [4] " 16 , 31008*mu [4] "15, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0]; 

> nops(pro20) ; 

> ccA [1 , 5] : =convert (pro20 , Vector) : 

> ccA[l,5] [55] ; 

pro20 := [0, 0, 0, 0, 0, 0, 1938 , 31008 /if , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0] 

65 
1 

> pro21 := [0,0, 0,0, 0,0, 320*mu [4] " 14 , 4480*mu [4] "13, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0]; 

> nops(pro21) ; 

> ccA[2,5] :=convert(pro21, Vector) : 

> ccA[2,5] [56] ; 

pro21 := [0, 0, 0, 0, 0, 0, 320 /ii"*, 4480 /if , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0] 

65 
1 

> pro22 := [0,0, 0,0, 0,0, 52*mu [4] " 12 , 624*mu [4] "11, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0]; 

> nops(pro22) ; 

> ccA [3 , 5] : =convert (pro22 , Vector) : 

> ccA[3,5] [57] ; 

pro22 := [0, 0, 0, 0, 0, 0, 52 /if , 624 /if , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0] 

65 
1 

> pro23 := [0,0, 0,0, 0,0, 8*mu [4] " 10 , 80*mu [4] "9, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0]; 

> nops(pro23) ; 

> ccA [4 , 5] : =convert (pro23 , Vector) : 

> ccA[4,5] [58] ; 

pro23 := [0, 0, 0, 0, 0, 0, 8 /if , 80 /t^, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0] 

65 
1 



> pro24:=[l/mu[l] "5,-5/mu[l] "6,0,0,0,0,mu[4] '8,8*mu[4] "7,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0]; 

> nops(pro24) ; 

> ccA [5 , 5] : =convert (pro24 , Vector) : 

> ccA[5,5] [59] ; 

pro24 ~ [/^J;^ -5 0' 0. 0' 0- M4> 8 o> 0, 0, 0. o> o> 0, 0, 0, o> 0> 0, 0, 0, o> 0> 0, 0, 0, o> 0> 0, 0, 0, 0> 0> 0: 0, o> 0> 0> 0: 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0] 

65 

1 

> pro25 := [0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0]; 

> nops(pro25) ; 

> ccA [6 , 5] : =convert (pro25 , Vector) : 

> ccA[6,5] [60] ; 

pro25 ~ [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,1,0, 0, 0, 0, 0] 
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Columns of yyj' {twenty Eve (10 x l)-vectors) 

> Mg:=(y,z,p,q)->(y,z+t*y"p/z"q) ; 

Mg := {y,z,p,q) ^y,z+ — 

> M: = (y,z)->Mg(y,z,p,q) ; 

M {y,z) Mg{y,z,p,q) 

> cccA:=unapply(MatrixMatrixMultiply(R_2,champ2(M)) , (p,q)) : 

Glueing 

We form the 25 vectors of type c of 1/(232) ■ 

> for q from 1 to 5 do 

> for p from to q+1 do 

> for j from 1 to 10 do 

> uu[p,q] [j] :=cA(p,q) [j] 

> od 

> od 

> od: 

> for q from 1 to 5 do 

> for p from to q+1 do 

> for j from 11 to 75 do 

> uu[p,q] [j] :=ccA[p,q] [j-10] 

> od 

> od 

> od: 

> for q from 1 to 5 do 

> for p from to q+1 do 

> for j from 76 to 85 do 

> uu[p,q] [j] :=cccA(p,q) [j-75] 

> od 

> od 

> od: 

> for q from 1 to 5 do 

> for p from to q+1 do 

> uu[p,q] :=convert(uu[p,q] ,list) 

> od 

> od: 

> q:=l: 

> for p from to q+1 do 

> VectDef _c [p+1] : =uu [p , q] 

> od: 



> q:=2: 

> for p from to q+1 do 

> VectDef _c [p+4] : =uu [p , q] 

> od: 

> q:=3: 

> for p from to q+1 do 

> VectDef _c [p+8] : =uu [p , q] 

> od: 

> q:=4: 

> for p from to q+1 do 

> VectDef _c [p+13] : =uu [p , q] 

> od: 

> q:=5: 

> for p from to q+1 do 

> VectDef _c [p+19] : =uu [p , q] 

> od: 

> p:='p': 

> q:='q': 



6.6.4 Vectors of type d (we compute the columns of the matrices yy^, ,Ad and ,yVJ^') 
Columns of (Eve (10 x l)-vectors) 

> Mg:=(y,z,p)->(y+t*y~(p+2)/z"p,z+t*y"(p+l)/z"(p-l)) ; 

Mg := [y,z,p) ?/H —rir 

zP zP 

> M:=(y,z)->Mg(y,z,p); 

M := {y,z) ^ Mg{y,z,p) 
Adp 

> dA:=unapply(MatrixMatrixMultiply(R_l,champl(M)) ,p) : 

Columns of yV^ to which we add 50 zeros then a (5 x 5)-identity block 
Adi 

> prol:= [0,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0]; 

> nops(prol) ; 

> ddA [1] : =convert (prol , Vector) : 

> ddA [1] [61]; 

prol := [0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,0] 
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Ad2 

> pro2 : = [0 , , mu [2] "2 , 2*mu [2], 0,0, 1,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0]; 

> nops(pro2) ; 

> ddA [2] : =convert (pro2 , Vector) : 

> ddA [2] [62] ; 

pro2 := [0,0,/^^, 2 ^2, 0,0, 1,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0] 

65 
1 

Ad3 

> pro3:=[0,0,0,0,mu[3] ~2,2*mu[3] ,2*mu[4] ,2,2*mu[5] ,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0]; 

> nops(pro3) ; 

> ddA [3] : =convert (pro3 , Vector) : 

> ddA [3] [63] ; 



pro3 [0, 0, 0, 0, ^^3, 2 /^a, 2 /i4, 2, 2 /ig, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0] 
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> pro4 : = [0 , , , , , , mu [4] " 2 , 2*mu [4], 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0]; 

> nops(pro4) ; 

> ddA [4] : =convert (pro4 , Vector) : 

> ddA [4] [64] ; 

pro4 := [0, 0, 0, 0, 0, 0, fil, 2 ^4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,1,0] 
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Ads 

> pro5 : = [0 , , , , , , , , mu [5] " 2 , 2*mu [5], 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1]; 

> nops(pro5) ; 

> ddA [5] : =convert (pro5 , Vector) : 

> ddA [5] [65] ; 

pro5 := [0,0,0,0,0,0,0,0,/i^, 2 ^5, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1] 
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Columns of (Eve (10 x l)-vectors) 

> Mg : =(y , z ,p) -> (y+t*y- (p+2) /z'p , z+t*y" (p+1) /z" (p-1) ) ; 

Mg := [y,z,p) ^'^+ ^TT 

zP zP ^ 

> M: = (y,z)->Mg(y,z,p) ; 

M {y,z) Mg{y,z,p) 
Adp 

> dddA:=unapply(MatrixMatrixMultiply(R_2,champ2(M)) ,p) : 
Glueing 

We form the 5 vectors of type d of V{1)2)- 

> for k from 1 to 5 do 

> for j from 1 to 10 do 

> VectDef_d[k] [j] :=dA(k) [j] 

> od 

> od: 

> for k from 1 to 5 do 

> for j from 11 to 75 do 

> VectDef_d[k] [j] :=ddA[k] [j-10] 

> od 

> od: 

> for k from 1 to 5 do 

> for j from 76 to 85 do 

> VectDef_d[k] [j] :=dddA(k) [j-75] 

> od 

> od: 

> for k from 1 to 5 do 

> VectDef_d[k] : =convert (VectDef _d [k] ,list) 

> od: 

> for k from 1 to 5 do 

> VectDef _d [k] : =convert (VectDef _d [k] , Vector) 

> od: 



6.6.5 Construction of the matrix '1^2 



We put together the 8 vectors of type a, the 25 vectors of type b, the 25 vectors of type c and the 5 vectors of 
type d. 

> for k from 1 to 8 do 

> V[k] :=VectDef_a[k] 

> od: 

> for k from 9 to 33 do 

> V [k] : =VectDef _b [k-8] 

> od: 

> for k from 34 to 58 do 

> V [k] : =VectDef _c [k-33] 

> od: 

> for k from 59 to 63 do 

> V [k] : =VectDef _d [k-58] 

> od: 

> for k from 1 to 63 do 

> WW[k] :=convert(V[k] ,list) 

> od: 

> MatVectDef :=Transpose(Matrix([seq(WW[k] ,k=l. .63)])) : 

We define the images of the 8 vectors of V^(Si) under the action of . 

> for k from 1 to 8 do 

> Image VectDef [k] : =MatrixMatrixMultiply(MM, VectDef [k] ) 

> od: 

We know how to compute these 8 vectors in the basis with the 65 vectors of V{D2) thanks to the relations (5.3) 
given in the article. We check that these relations are satisfied (these 8 instructions give 85 zeros). 

> for k from 1 to 85 do 

> simplifyClmageVectDef [1] [k]-WW[l] [k]-3*WW[62] [k]+4*WW[15] [k]+3*WW[36] [k] ) ; 

> od: 

> for k from 1 to 85 do 

> SimplifyClmageVectDef [2] [k]-WW[2] [k]+3*WW[60] [k] ) ; 

> od: 

> for k from 1 to 85 do 

> SimplifyClmageVectDef [3] [k]-WW[3] [k]+3*WW[59] [k] ) ; 

> od: 

> for k from 1 to 85 do 

> SimplifyClmageVectDef [4] [k]-WW[4] [k]-2*WW[20] [k]-WW[40] [k]+2*WW[63] [k] ) ; 

> od: 

> for k from 1 to 85 do 

> SimplifyClmageVectDef [5] [k]-WW[5] [k]-2*WW[61] [k] ) ; 

> od: 

> for k from 1 to 85 do 

> SimplifyClmageVectDef [6] [k]-WW[6] [k]-2*WW[60] [k] ) ; 

> od: 

> for k from 1 to 85 do 

> SimplifyClmageVectDef [7] [k]-WW[7] [k] ) ; 

> od: 



6.7 Computation of the action of ip^, on H^(X, TX) 
Construction of the supplementary (fi of V{'Di) of dimension 22 in W{Di). 

We form a (30 x 30)-matrix NNi whose the first 22 columns are now zero and the last 8 are the basis vectors of 
V{Di) computed previously. 

> for k from 1 to 30 do 

> for j from 23 to 30 do 

> Nl [k] [ j ] : =VectDef [ j -22] [k] 

> od 

> od: 

> for k from 1 to 30 do 

> for j from 1 to 22 do 

> Nl[k][j]:=0 

> od 

> od: 

> for k from 1 to 30 do 

> Nl [k] :=convertCNl [k] ,list) 

> od: 

> NNl:=MatrixC[seqCNl[k] ,k=l. .30)]): 



We fill the first 22 columns of NNi by some vectors chosen among the algebraic basis oiW{'Di); they will define 
a supplementary (fi of y(J)i). 

> NN1[1,1]:=1: 

> NN1[2,2]:=1: 

> NN1[3,3]:=1: 

> NN1[4,4]:=1: 

> NN1[5,5]:=1: 

> NN1[6,6]:=1: 

> NN1[7,7]:=1: 

> NN1[8,8]:=1: 

> NN1[9,9]:=1: 

> NN1[10,10] :=1: 

> NN1[11,11] :=1: 

> NN1[12,12] :=1: 

> NN1[13,13] :=1: 

> NN1[14,14] :=1: 

> NN1[15,15] :=1: 

> NN1[16,16] :=1: 

> NN1[17,17] :=1: 

> NN1[21,18] :=1: 

> NN1[22,19] :=1: 

> NN1[23,20] :=1: 

> NN1[24,21] :=1: 

> NN1[25,22] :=1: 

We compute the determinant of NNi, that is nonzero for generic parameters fii] this shows that S'l is a supple- 
mentary to V{^i). 

> Determinant (NNI ) ; 

12 ^4/i5(— 6 fi4fi5 — 4 fi4fil + 9 fi^a + 9 /^a + 45 ^5 + 45 fil — 36 ^14) 

We extract the first 22 columns of NNi in order to form the matrix left, that is the injection matrix of S'l into 

> for k from 1 to 30 do 

> for j from 1 to 22 do 

> coeff 1 [k] [j] :=NN1 [k, j] : 

> od: 

> od: 

> for k from 1 to 30 do 

> coeff 1 [k] : =convert (coeff 1 [k] , list) 

> od: 

> left:=Matrix( [seq(coeff 1 [k] ,k=l . .30)] ) : 

Definition of the supplementary (02 of V{D2) in W{'I>2) of dimension 22 which is the image of #1 by the inclusion 
ofW{^i) into W{D2). 

We form a (85 x 85)-matrix NN2 whose the first 22 columns are now zero and the last 63 are some basis vectors 
of V{^2) computed previously. 

> for k from 1 to 85 do 

> for j from 23 to 85 do 

> N2[k] [j] :=WW[j-22] [k] 

> od 

> od: 

> for k from 1 to 85 do 

> for j from 1 to 22 do 

> N2[k][j]:=0 

> od 

> od: 

> for k from 1 to 85 do 

> N2[k] :=convert(N2[k] ,list) 

> od: 

> NN2:=Matrix( [seq(N2[k] ,k=l. .85)]): 

We fill the first 22 columns of NN2 by the basis vectors of (02 . 

> NN2[1,1]:=1: 

> NN2[2,2]:=1: 

> NN2[3,3]:=1: 

> NN2[4,4]:=1: 

> NN2[5,5]:=1: 



> NN2[6,6]:=1: 

> NN2[7,7]:=1: 

> NN2[8,8]:=1: 

> NN2[9,9]:=1: 

> NN2[10,10] :=1: 

> NN2[11,11] :=1: 

> NN2[12,12] :=1: 

> NN2[13,13] :=1: 

> NN2[14,14] :=1: 

> NN2[15,15] :=1: 

> NN2[16,16] :=1: 

> NN2[17,17] :=1: 

> NN2[76,18] :=1: 

> NN2[77,19] :=1: 

> NN2[78,20] :=1: 

> NN2[79,21] :=1: 

> NN2[80,22] :=1: 

The matrix right is the projection matrix of W{D2) onto £2. 

> NNN2 : =MatrixInverse (NN2) : 

> for k from 1 to 22 do 

> for j from 1 to 85 do 

> coeff2[k] [j] :=NNN2[k, j] : 

> od: 

> od: 

> for k from 1 to 22 do 

> coef f 2 [k] : =convert (coef f 2 [k] , list) 

> od: 

> right:=Matrix([seq(coeff2[k] ,k=l. .22)]) : 

The matrix OMG is the matrix of acting on H^(X, TX). 

> OMG : =MatrixMatrixMultiply (right , MatrixMatrixMultiply (MM , left) ) 

> OMG; 

' 22 X 22 ' (Matrix) 
'Data Type: ' anything 
' Storage : ' rectangular 
' Order : ' Fortran _order 
We assign the parameters et /i^ otherwise the computation is too heavy. 

> lambda [1] : =10: 

> lambda [2] : =-6: 

> lambda [3] : =5 : 

> lambda [4] : =2 : 

> lambda [5] : =3 : 

> mu[l]:=3: 

> mu[2]:=-2: 

> mu [3] : =4 : 

> mu [4] : =3 : 

> mu[5]:=29: 

Computation of the eigenvalues of OMG. 

> f actor (CharacteristicPolynomial (OMG, x)) ; 

(a;2 + 3x + l){x^ + 18x + l){x^ - 7a; + l){x^ + x+ l){x - lf{x + lf{x^ ~ x 
We look at the size of the Jordan blocks for the multiple eigenvalues. 

> Id : =Matrix (22 , 22 , shape=identity) : 

> Rank(OMG+ld) ; 

> Rank((0MG+ld)"2) ; 

> Rank((0MG+ld)"3) ; 

19 
18 
18 



> Rank(OMG-Id) ; 

> Rank((0MG-Id)'2) ; 

20 
20 

> Rank(0MG-(l+I*sqrt(3))/2*Id) ; 

> Rank((0MG-(l+I*sqrt(3))/2*Id"2)) ; 

19 
19 

It takes 13 minutes to execute the file. 
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